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ABSTRACT 

We present a new deep 21-cni survey of the Andromeda galaxy, based on high resolution observations 
performed with the Synthesis Telescope and the 26-m antenna at DRAG. The Hi distribution and 
kinematics of the disc are analyzed and basic dynamical properties are given. The rotation curve is 
measured out to 38 kpc, showing a nuclear peak at 340 km s~^, a dip at 202 km s~^ around 4 kpc, 
two distinct flat parts at 264 km s~'^ and 230 km s~^ and an increase to 275 km in the outermost 
regions. Except for the innermost regions, the axisymmetry of the gas rotation is very good. A very 
strong warp of the Hi disc is evidenced. The central regions appear less inclined than the average disc 
inclination of 74°, while the outer regions appear more inclined. Mass distribution models by ACDM 
Navarro-Frenk- White, Einasto or pseudo-isothermal dark matter halos with baryonic components are 
presented. They fail to reproduce the exact shape of the rotation curve. No significant differences are 
measured between the various shapes of halo. The dynamical mass of M31 enclosed within a radius 
of 38 kpc is (4.7 ±0.5) x 10^^ A^o- The dark matter component is almost 4 times more massive than 
the baryonic mass inside this radius. A total mass of 1.0 x 10^^ A4q is derived inside the virial radius. 
New Hi structures are discovered in the datacube, like the detection of up to five Hi components 
per spectrum, which is very rarely seen in other galaxies. The most remarkable new Hi structures 
are thin Hi spurs and an external arm in the disc outskirts. A relationship between these spurs and 
outer stellar clumps is evidenced. The external arm is 32 kpc long, lies on the far side of the galaxy 
and has no obvious counterpart on the other side of the galaxy. Its kinematics clearly differs from 
the outer adjacent disc. Both these Hi perturbations could result from tidal interactions with galaxy 
companions. 

Subject headings: galaxies: ISM - galaxies: fundamental parameter (mass) - galaxies: individual 
(M31, NGC 224) - galaxies: kinematics and dynamics - galaxies: structure - Local 
Group 



1. INTRODUCTION 

In the last decade, several efforts have led to a better 
understanding of the building history of the disc of the 
Andromeda galaxy (hereafter M31, Table 1), the most 
nearby massive spiral galaxy to the Milky Way (MW). 

One of the most exciting result is probably the discov- 
ery of an extended stellar halo around M31 (Ibata et al. 
2001). More and more studies are completing this new 
view of the stellar distribution of M31 and are reveal- 
ing little by little the complex chemical, kinematical and 
morphological nature of its halo (Ferguson et al. 2002; 
Reitzel & Guhathakurta 2002; Ibata et al. 2004, 2005; 
Brown et al. 2006; Guhathakurta et al. 2006; Kalirai et 
al. 2006a,b; Martin et al. 2006; Faria et al. 2007; Ibata et 
al. 2007; Majewski et al. 2007; Brown et al. 2008; Chap- 
man et al. 2008). Basically, these authors find that the 
outskirts of the stellar disc are very clumpy and irregular. 
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and the extended halo is a mix of many metal poor and 
rich substructures like the south-eastern Giant Stream 
(GS) and other seemingly shorter streamlike extensions, 
low luminosity dwarf satellites and globular clusters, as 
well as an old metal poor underlying (primordial?) halo 
which extends up to at least 150 kpc (in projection). 
These stellar tidal features and companion galaxies are 
the probable main imprints of the hierarchical growth 
of the M31 stellar disc and halo, like the ones encoun- 
tered in numerical models of dark matter evolution in the 
framework of the Cold Dark Matter (CDM) paradigm 
(e.g. Frenk et al. 1988; Klypin et al. 1999; Springel et al. 

2005) . Numerical models trying to explain the presence 
of the GS tidal debris around M31 propose a tidal inter- 
action with a ^ 10^ A4q galaxy progenitor (Fardal et al. 

2006) . 

The gas content in the halo of M31 has also been 
probed recently from deep Hi measurements, allowing 
the detection of a population of discrete low luminos- 
ity clouds, with masses in the range of 10^ to 10^ Mq 
(Thilker et al. 2004; Westmeier, Braun & Thilker 2005). 
The authors suggest that those clouds could be the 
analogs of the high-velocity Hi clouds seen around the 
MW (Wakker et al. 1999). They also suggest a manyfold 
origin of that gas (residual from galaxy merger or inter- 
action, cooling flow in the Local Group, etc.), which 
may only contribute to ^ 1% of the total Hi disc mass of 
M31. Hi gas clouds from the intergalactic medium could 
thus also be the building blocks of M31. 

At the same time, the advent of the new generation of 
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X-ray, ultraviolet and infrared space observatories allows 
to probe the structure of the interstellar medium of M31 
with unprecedented details (Williams et al. 2004; Thilker 
et al. 2005; Barmby et al. 2006; Gordon et al. 2006; 
Stiehle et al. 2008). For instance, the mid- infrared im- 
ages reveal a dominant ring-like distribution for the dust 
component. Its perturbed morphology seems to have 
been shaped by the passage of the companion Messier 
32 through the M31 disc (Thilker et al. 2005; Block et 
al. 2006). A high resolution kinematical survey of the 
molecular gas contained in the ring and spiral structures 
has also been presented (Nietcn et al. 2006). 

A large part of the Hi studies that have been done on 
M31 go back up to ~30 years ago (Guibert 1973; Emer- 
son 1974; Newton & Emerson 1977; Cram, Roberts & 
Whitehurst 1980; Unwin 1980a,b; Brinks & Shane 1984). 
These data were acquired at low spectral and/or angular 
resolution data and did not necessarily cover the whole 
disc of the galaxy, mainly for sensitivity reasons. An- 
other more recent Hi study of M31 has been presented in 
Braun (1990) from VLA observations. Perhaps the most 
important result from all these studies is that the Hi disc 
of M31 exhibits a warp whose effects are remarkable in 
the datacube by the presence of many spectral Hi peaks 
along different line-of-sights. Another interesting result 
by Braun (1991) is the modeling of a rotation curve that 
appears to decline as a fimction of the galactic distance. 
A study of the mass distribution using this rotation curve 
only requires stellar bulge plus disc components, with no 
need of any dark matter halo. If it is really the case, M31 
would be very different from every other spiral galaxies 
which are known to exhibit a flat rotation curve (or even 
increasing) at large galactocentric distances and thus to 
contain a massive hidden mass, unless the law of Gravity 
is modified in these acceleration regimes (Milgrom 1983, 
2008). 

This result was the main motivation to get single dish 
observations along the photometric major axis for the 
approaching disc half. Indeed, radial velocities of the Hi 
gas in the receding side of M31 are contaminated by Hi 
gas in the MW. We concluded that the outer Hi rotation 
curve can not be decreasing but seems remarkably flat at 
the largest radii (Carignan et al. 2006). However, since it 
was not possible to model the warp properly with those 
single dish data, it was decided in 2005 to get full 2D ve- 
locity information from wide-field synthesis observations 
at DRAG. 

In this article, we present preliminary morphological, 
kinematical and dynamical results obtained from the 
analysis of the DRAG Hi observations. The direct ob- 
jectives of this article are to present the most extended 
Hi distribution of M31, to derive an accurate Hi rotation 
curve in order to verify whether the rotation velocities 
really decrease at large galactocentric distances or really 
remain constant, as well as to derive its basic dynamical 
parameters. 

It is worth mentioning here the very recent, wide- field 
and high angular Hi imaging presented in Braun et al. 
(2009) with the help of the WSRT and GBT telescopes. 
Future results from that very deep dataset will surely 
serve as independent comparison with those we present 
in this article and in other forthcoming papers from this 
series. 

The article is organized as follows. Section 2 describes 



TABLE 1 
Basic optical parameters of M31. 
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the DRAG observations and the basic data reduction 

steps. Section 3 describes the line fitting procedure as 
well as how the Galactic Hi was subtracted from the 
datacube. Section 4 describes the Hi content and dis- 
tribution while section 5 concentrates on the kinematics 
and on the calculation of the Hi rotation curve. A com- 
parison with results from different studies is also done in 
this section. An analysis of the perturbed outer regions 
of the disc is done in section 7. Finally, a study of mass 
distribution models is done in section 8. Concluding re- 
marks are given in section 9. 

All velocities are given in the heliocentric rest frame 
and a Bubble's constant of Hq = 73 km Mpc~^ is 
chosen throughout the article (Spergel et al. 2007). Since 
the symbol R refers to as the galactocentric radius and 
most of the magnitudes and luminosities are given in the 
photometric i?— band, no subscripts are attached to mag- 
nitudes or luminosities for clarity reasons, except where 
explicitely mentioned. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. HI emission line observations 

A total of five fields towards M31 were observed in 
the 1420 MHz continuum and 21 cm line with the Syn- 
thesis Telescope (ST) at the Dominion Radio Astro- 
physical Gbservatory (DRAG) between September and 
December 2005. This interferometer consists of seven 
9 metres diameter antennae along an East- West baseline 
617.1 m long. The primary beam of each element is 107'. 2 
(FWHM) , and structures down to the resolution limit of 
58" X 58"/ sin((5) in the 1420 MHz line are resolved within 
this beamwidth (a Gaussian taper to the u, v data is ap- 
plied, broadening somewhat the synthesized beamwidth 
up from the nominal 49" x 49"/ sin(^) resolution in the 
continuum). Further instrumental details on the DRAG 
ST are found in Landccker et al. (2000). 

Tables 2 and 3 list the observational parameters and 
field-of-view centres for the five individual aperture syn- 
thesis observations. Fields were chosen such that a spac- 
ing of A =77' between centres was observed, giving 
nearly equal sensitivity over the whole area surveyed. 
The exposure times per field was 144 hours. 
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TABLE 2 

Observational parameters of the DRAO aperture 

SYNTHESIS observations OF MESSIER 31. 



Parameter 


Value 


Observation dates 


2005 September-December 


Total length of observation 


144x5 hrs 


Velocity centre of band 


-300 km s-i 


Total bandwidth 


4 MHz (843 km s-i) 


Number of velocity channels 


256 


Frequency sampling 


15.6 kHz 


Velocity resolution 


5.3 km 


Number of spatial pixels" 


1024 


Pixel angular sizc"^ 


22" 



^Values given for the full resolution initial datacube. 



TABLE 3 

DRAO Synthesis Telescope observational set-up and noise. 



Obs. 


Field centre 


Beam parameters 

^maj ( ) •'^ ^min ( ) 

and orientation (°) 


Noise at field 


date 

(2005) 


coordinates 
(J2000, h:m:s,°:':") 


centre 
(K) 


Oct. 05 


0:52:19.7, 43:15:00 


1.46x0.98, -90.2 


1.12 


Sep. 09 


0:47:32.0, 42:19:00 


1.50x0.98, -90.2 


0.98 


Sep. 09 


0:42:44.3, 41:16:09 


1.53x0.98, -90.0 


0.83 


Nov. 04 


0:37:56.6, 40:08:00 


1.58x0.97, -89.8 


1.04 


Dec. 10 


0:33:18.9, 39:00:00 


1.57x0.97, -87.2 


1.06 



Summary of 1420 MHz line + continuum observations centered 
on and surrounding Messier 31. Noise is indicated for the full 
resolution initial datacube. 

Hi emission line images arc made in each of 256 chan- 
nels across a 4 MHz bandwidth centered onv — —300 km 
s~^. The spectrometer gives a resolution of 5.27 km s~^, 
and each channel samples a width of Av = 3.3 km s~^. 
The Hi line datacube spans 843 km s~^. The theoretical 
noise at the pointing centre is 1.75 sin((5) K, correspond- 
ing closely to the measured values in Table 3 (for channels 
free from Hi emission). 

To depict structures accurately in the radio continuum 
and Hi line images (especially those of angular size > 56', 
missed by the interferometer due to its shortest base- 
line limit of 12.9 m), low spatial frequency information is 
routinely added to each ST map. These "short-spacing" 
data are obtained from Hi line observations made with 
the DRAO 26-m paraboloid (Higgs & Tapping 2000). 
The flux is corrected for stray radiation entering through 
the sidelobes, and the continuum is subtracted from all 
channels before the integration with the ST observations. 
The spatial resolution of these data is 37' x35.'3; all spec- 
trometer settings were the same as for the ST data. 

2.2. Data Reduction and Mosaicing 

Prior to mosaicing our five fields together, we per- 
form some standard data processing steps (processing 
methods arc similar to those developed for the Cana- 
dian Galactic Plane Survey CGPS, described in Taylor 
et al. 2003). First, a continuum baseline level is removed 
from each datacube. Average emission line-free channels 
at both the low and high velocity ends of the cube are 
made, and a linear interpolation (determined from these 
two continuum maps) is then subtracted from each chan- 
nel map in the datacube. To calibrate the flux scale of 
each cube, we compare point sources in the average of 
the two continuum end-channels to those in the 30 MHz 
continuum band map of each field (these maps are first 
CLEANed around the strongest sources). This can be 



done since continuum maps are flux calibrated against 
several strong sources that arc routinely observed by the 
ST (e.g. 3C 48, 3C 286; see Table 1 of Taylor et al. 2003, 
for sources and fluxes). An error- weighted mean of the 
flux ratio 51420 coih/'Shi cont for all point-sources above a 
cutoff level (20 mJy beam~^) is obtained; this sample is 
further trimmed by rejecting sources ±0.5 away from the 
mean. Typically, 20-30 sources remained in the sample 
for each field. The uncertainty in the flux calibrated this 
way is ~5%. 

The processed fields are combined into a 1024 square 
pixels mosaic, centered on optical coordinates RA = 
00M2'"44.3", DEC = +41° 16' 09" (J2000). The pixel 
size is 21'.'875. The central 93' radius of the primary 
beam constitutes the final width of the five individual 
fields, which were mosaiced together. The final mosaic 
spans 515' oriented along M31's disc. The final noise val- 
ues for the mosaic (measured) are AT;, ~0.85 K in the 
individual field centres, and ATf, ^0.95 K in the overlap- 
ping regions. 

The final step consists in merging interferometer data 
with single-dish (26-m telescope) observations. For sim- 
plicity, this merging is performed in the map plane (as 
opposed to the u, v visibilities plane, as is done for the 
CGPS; see Taylor et al. 2003) using the procedure de- 
scribed below. The interferometer mosaic is first con- 
volved to the resolution of the single-dish mosaic; this is 
then subtracted from the 26-m mosaic. The resulting dif- 
ference map should contain Hi structures invisible to the 
interferometer, so it is added back to the full-resolution 
ST mosaic. This procedure gives equal weight to low and 
high spatial frequency structures, and while it is some- 
what different from the method used in the CGPS (where 
a tapering function weights the overlapping spatial fre- 
quencies of both maps), the results of both approaches 
are very similar for Hi line images. 

A spatial binning of 2x2 pixels is finally applied to 
the datacube in order to decrease noise in the spectra. 
The datacube has a final dimension of 512x512 pixels, 
with a pixel size of 43'.'75. This angular sampling is well 
sufficient for the kinematical and dynamical purposes. It 
corresponds to 167 pc at the distance of 785 kpc. The 
DRAO synthesized beam size of ^ 60" x 90" samples a 
linear scale of ^ 230 x 340 pc in the Hi disc of Messier 31. 
It is as resolved as a typical Hi mapping of 10 Mpc distant 
discs observed with a high resolution beam resolution 
of e.g. 5". The present Hi observations can thus be 
considered as high resolution ones and for this reason no 
correction for beam-smearing is needed. 

3. Hi datacube ANALYSIS 

3.1. Previous works 

Hi spectra of Messier 31 are known to exhibit multiple 
peaks in emission (e.g. Newton & Emerson 1977; Cram, 
Roberts & Whitehurst 1980; Bajaja & Shane 1982). It is 
also the case in molecular gas observations (Dame et al. 
1993; Loinard et al. 1999; Nieten et al. 2006). For the Hi 
data, this is particulary well illustrated along position- 
velocity (PV) diagrams made parallel to the major axis in 
Cram, Roberts & Whitehurst (1980) or Brinks & Burton 
(1984), where a second velocity peak is seen in addition 
to the main (brighter) peak. Both these lines draw a 
tilted figure-eight shape in which the main peak traces 
the usual rotating pattern of the disc (steep velocity rise 




in the inner regions sometimes accompanied by a flatter 
velocity part in the external regions) and the second peak 
shows a linear (shallower) velocity rise all along the slices. 

With the noticeable exceptions of Brinks & Shane 
(1984), Brinks & Bm-ton (1984) or Braun (1991), the 
previous Hi studies of M31 analyzed their data with a 
single emission line approximation in order to derive an 
integrated emission map, a velocity field or a rotation 
curve. For instance, Newton & Emerson (1977) reported 
that when two peaks are detected, the radial velocity of 
the spectrum is chosen at the barycentre of the lines. 
Though it seems a reasonable hypothesis in regard to 
the low spectral resolution of all these old observations, 
it surely provides a biased velocity distribution for M31. 
Integrated fluxes and velocities are uncorrectly estimated 
when lines are blended. 

A simple explanation for the origin of two components 
comes from the fact that the Hi disc of Messier 31 is 
warped at large galactocentric radius. Because the disc 
is highly inclined, the projection effects enable us to cross 
two times the disc along the line-of-sight (Hi gas is sup- 
posed to be optically thin). The line-of-sight velocities 
are thus composed of a main Hi component lying in the 
internal disc region in addition to a secondary compo- 
nent lying in the external warped region of the disc but 
seen projected at small radii. The presence of the warp 
were the reasons that led Brinks & Burton (1984) to an- 
alyze the WSRT data differently from the other studies 
by a warped and flaring Hi layer model. They proposed 
two velocity fields, one for the disc and another one for 
a separate, external warped structure. From this, they 
deduced that about 39% of the total Hi mass of the disc 
could reside beyond 18 kpc, in its warped region. 

3.2. Current analysis 
3.2.1. Evidence for multiple HI spectral components 

A provisional study of the DRAO datacube has been 
presented in Carignan, Chemin & Foster (2007) where 
moments maps were derived by fitting a single gaussian 
to each spectrum of the datacube. We extend here this 
simplistic view of our data to a more detailed analy- 
sis where multiple Hi components are now fitted to the 
spectra. The high-sensitivity observations indeed reveal 
a new more complex view of the neutral gas in M31 
than what has been shown in past Hi studies. Figure 1 
displays several Hi spectra picked at different positions 
through the datacube. The Hi spectra exhibit almost 
all the time more than one component, sometimes up to 



5 peaks. This result confirms the basic view shown in 
e.g. Brinks & Burton (1984) where two Hi components 
are shown (and sometimes more) , but largely extends the 
detection of new multiple spectral components in the dat- 
acube. One clearly wants to emphasize here that a single 
or two peak analysis cannot apply straightforwardly to 
the current data. This would provide biased intensity, 
velocity and velocity dispersion maps. 

3.2.2. Contamination by the Milky Way 

A usual problem with observations of M31 is the con- 
tamination from the MW. M31 is so massive that its re- 
ceding half can reach radial velocities that coincide with 
those of the Galaxy. However most of Galactic Hi re- 
mains relatively easy to detect and to model because its 
lines can be observed almost everywhere in the field of 
view (Fig. 2). We detect two major Galactic lines. A 
first bright component is found at — 4 ± 2 km s~^ for 
a velocity dispersion of 6 ± 3 km s^^ as deduced from 
66000 spectra. This line does not contaminate the emis- 
sion from Messier 31 due to its too high radial velocity. 
Another line is found at —40 ±11 km s^^ for a veloc- 
ity dispersion of 18 ± 10 km s^^, as deduced from 19500 
spectra free from M31 Hi emission. This is the line that 
mostly contaminates the NE half of M31. Another mi- 
nor Milky Way line is observed around —20 km s~^ (see 
Fig. 2). This emission line does not contaminate the sig- 
nal from M31 because of its high radial velocity. Other 
very faint lines only revealed by very high contrast imag- 
ing are found around —70 km s~^ and —90 km s~^. Their 
contribution to the whole emission is very negligible. 

The method used to subtract the MW emission from 
the datacube is similar to the one employed by Braun 
et al. (2009). For each channel map, a mask containing 
pixels with Hi emission of M31 is first created. The mask 
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Fig. 2. — Contamination of the data by Galactic Hi. The slice 
through the datacube shows 3 Hi components from the Milky Way 
at ~ -4 km s'^ (1), ~ -20 km s"! (2) and ~ -40 km s"! (3). 
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the profile. 

covers a small fraction of the field-of-view and allows to 
blank the Hi emission of M31. A two-dimensional model 
of the Galactic emission is then generated by fitting a 
S"^*^ order 2D-polynomial to the blanked channel map. 
The 2D model is then subtracted from the initial channel 
map. This method allows to remove most of the galactic 
Hi. Residual MW contamination is removed by hand 
during the cleaning of the multiple velocity fields (see 
§3.2.3). We tried another approach to remove MW Hi 
by simultaneously fitting gaussian lines in addition to Hi 
from M31. However that subtraction method had to be 
rejected because too many negative residuals are created 
in the cube, introducing artificial local minima in the 
integrated emission map. 

As claimed in Braun et al. (2009) it is very difficult to 
define whether small scale Hi structures observed around 
the disc of M31 are bound to Andromeda or simply 
Galactic cloudlets. For simplicity reasons, we decide to 
ignore the gas emission out of the main disc in the fol- 
lowing analysis. 

There is no clear mention in older Hi studies of M3f 
of Galactic emission around ~ —40 km s~^. It is thus 
possible that measurements of the Hi flux and velocity 
dispersions are overstimated for the receding half of the 
disc in those works. However their velocity fields remain 
very marginally contaminated. 

3.2.3. Emission line detection and fitting 

The detection of emission lines is done by keeping all 
channels above a 3(t level, where a is the intensity dis- 
persion derived from half of the spectrum. This disper- 
sion measurement is justified by the fact that most of the 
channels are free from Hi emission. The method to locate 
a peak in the remaining channels is a usual slope compu- 
tation technique. An emission line is detected when the 
slope sign does not change along at least two consecutive 
channels, as measured both to the left and right sides 
of a peak. Once detected, velocities and amplitudes of 
these peaks serve as initial guesses for the minimization 
step where gaussian profiles are fitted to the spectrum. 
A "continuum" value is fitted as an additional parame- 
ter as well. It is very close to zero because a baseline 



has been subtracted in each pixel during the data re- 
duction (see §2.2). Boundary constraints are applied to 
the parameters (amplitude, velocity centroid and disper- 
sion). A velocity centroid has to be fitted within the 
observed spectral range, a velocity dispersion has to be 
greater than one channel width and lower than an artifi- 
cial (unexpected) large value of 130 km s~^ and a peak 
amplitude has to be greater than zero and lower than 
the highest intensity observed in a spectrum. Results of 
spectral fittings of 1 to 5 Hi components are displayed in 
Fig. I. 

A final cleaning step is then applied to the data. All 
fitted lines whose amplitude is larger than a threshold 
of 3(7 were kept. All lines whose velocity dispersion is 
greater than 60 km s~^ are discarded because such high 
values turn out to be very rare and, above all, unrealistic. 
Obvious residuals of MW emission are also removed from 
the different maps. As a result of this filtering procedure, 
among the 25350 available spectra in the Hi datacube 
and in which at least one component is detected with a 
signal-to-noise ratio of at least 3, two peaks are detected 
in ~66% of them, three peaks in ~35% of them, four 
peaks in ~14% of them and five peaks in ~3% of them. 
A large number of Hi components (> 3 peaks) is therefore 
rare in the datacube. 

4. HI DISTRIBUTION AND KINEMATICS 

4.1. Integrated properties 

The flux in each individual channel was summed to give 
the global Hi profile of Figure 3. The profile is relatively 
regular. The two peaks are almost perfectly symmetric. 
An asymmetry is nonetheless observed between them. 
There is more gas in the receding side of the galaxy than 
in its approaching part. An intensity- weighted systemic 
velocity of —299 ± 3 km s~^ and a midpoint heliocentric 
radial velocity oi v = —297 ± 3 km s~^ are derived. The 
measured profile widths at 20% and 50% levels are W20 = 
533 ± 3 km s'^ and W50 = 509 ± 3 km 8"^ This can 
be compared to the RC3 (de Vaucouleurs et al. 1991) 




45™ 40° 

Right Ascension 

Fig. 4. — Total Hi distribution of Messier 31. The small ellipse in 
the bottom-left corner displays the angular size of the synthesized 
beam (~ 1' X 1.5'). 
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Fig. 5. — 3D view of the Hi datacube of M31. The middle panel 
is the position- velocity diagram of the whole datacube projected on 
the photometric major axis, which is aligned along the horizontal 
axis. The bottom panel is the total integrated Hi image. The top 
panel is the mirror image of the bottom panel. Solid (dashed) lines 
are for Hi emission from the front (respectively far) side of the 
galaxy. Labels Nl, N2, SE show the two northern Hi spurs and 
the south-western extension, EA the external arm, Rl and R2 the 
ring-like structures in the centre of the disc (see text for details, 
§4.2). 

values oiv = -300 ± 4 km s-\ W20 = 536 ± 7 km s"! 
and W50 = 510±7 km s~^, which show a good agreement 
with ours within the errors. 

The integrated flux of 29221.4 Jy km s~^ implies a to- 
tal Hi mass of 4.23 x 10^ Mq. This mass is larger by 
15% than the one measured in Brinks & Burton (1984), 
is comparable with the estimate by Cram, Roberts & 
Whitehurst (1980) and is smaller by 25% than the one 
derived in Braun et al. (2009) (before opacity correc- 
tions). The difference between all mass estimates is likely 
due to the various sensitivities of the observations and to 
the accuracy of the subtraction of the contamining MW 
emission. 

4.2. The 3D structure of M31 

Figure 4 displays the total Hi integrated emission of 
M31 and a 3D view of the datacube is shown in Fig. 5. 
The lowest Hi column density of the DRAO observations 
at a 3(7 peak detection level is ^ 1.7 x 10^^ cm~^ and 
the highest Hi column density is 5 x 10^^ cm^^. The 
maximum disc extent is R ~2.6° (i? = 36.4 kpc), as 
derived along the photometric major axis. It corresponds 
to 1.67 i?25 and 6.48 with Rd being the stellar disc 
i?— band scale-length derived in §8.2.1. 

The high resolution Hi map reveals a disc with very lit- 
tle gas in its central regions, a feature usually observed 
in early- type discs. Two ring-like structures are observed 
around R - 10' (2.5 kpc) and R - 20' -25' (4.6-5.7 kpc). 
We refer to them as Rl and R2 in Fig. 5. At such high 
inclination it is difiicult to firmly claim whether they are 



real rings, like e.g. those created by gas accumulation 
at the location of inner or ultra-harmonic Lindblad res- 
onances, or tightly wound spiral arms. They are coin- 
cident with dusty ring-like structures observed in NIR 
images from Spitzer/IRAC data (Barmby et al. 2006) as 
well as with molecular gas ring-like structures (Nieten 
et al. 2006). Another wider bright ring- like structure 
is observed for 40' < i? < 80' (9.1-18.3 kpc). Though 
it is often referred to as the "lO-kpc ring" of M31, its 
morphology is more complex than a regular ring because 
holes are observed in it, as seen e.g. in the approaching 
side of the disc. Long spiral arms are then clearly seen 
at large radii. The most prominent spiral arm of M31 is 
observed in the south-western half of the disc. 

A faint external spiral arm (label EA of Fig. 5) is 
discovered on the edge of the receding half of the disc 
(southern part of the NE quadrant of Fig. 4). It is 
connected to the long spiral arm which arises from the 
SW and its apparent end is clumpy (see region around 
a2ooo = 00'^49'"41.80^ (52000 = +42'^12'54"). Other new 
structures that were not seen in old Hi data are the 
two disc extremities to the SW and NE (labels "SE" 
for south-western extension and "N1"/"N2" for northern 
spurs in Fig. 5). Here thin and faint gaseous extensions 
are observed. Their gas distribution and kinematics are 
discussed in more details in §7. 

The top and bottom panels of Fig. 5 is the right 
ascension-declination view of the datacube while the 
middle panel is a rotation of this later by 90° with respect 
to the horizontal axis (disc photometric major axis). It 
is thus a position velocity plot of the whole datacube 
projected onto the major axis. Vertical lines are drawn 
to guide the eye in order to link the kinematics to the 
spatial distribution of the Hi gas. The main features of 
the diagram are : 

• All spiral- and ring-like structures do not cross at 
the same location in the centre of the diagram 
{v — —300 km s^^). This is probably caused by 
a lopsided nature of the disc. 

• Except for the external arm (EA), all spiral- and 
ring-like structures have a symmetric counterpart 
with respect to the galaxy centre. This symmetry 
is not regular in velocity amplitude within ~ 50' 
around the centre, as seen for instance with the 
velocity peaks of the ring- like structure R2. This 
is also caused by the disc lopsidedness, as well as 
other probable noncircular gas motions in the cen- 
tral parts. 

• A steep velocity gradient is observed in the inner- 
most ring- like structure (Rl). 

• The northern spur N2 appears as a kinematical ex- 
tension of the external arm in the velocity space. 

The neutral gas distribution of the DRAO observa- 
tions is consistent with Hi images from many other works 
(Cram, Roberts & Whitehurst 1980; Roberts & White- 
hurst 1975; Bajaja & Shane 1982; Emerson 1974; Brinks 
& Shane 1984). The agreement is better with the deep 
WSRT observations of Braun et al. (2009) because of 
the comparable sensitivity with the DRAO observations. 
Their Hi image do show the external arm as well as the 
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north-eastern spurs and south-western extension. It thus 
leaves no doubt about their presence in the disc outskirts 
of M31. 

The gas distribution for all detected Hi components is 
shown in Figure 6. In each pixel, the components are 
displayed by decreasing integrated emission from left to 
right. A brief discussion on the origin(s) of all these Hi 
peaks is proposed in §6 

5. KINEMATIC AL ANALYSIS 

In Brinks & Burton (1984), a 3D model of WSRT ob- 
servations was presented in order to reproduce the com- 
plex structure of the datacube. They proposed a flared 
and warped disc model for M31 that can explain the 
presence of double Hi peaks observed in most of pixels. 
Not all of the dozen free parameters of their model are di- 
rectly fitted to the datacube. Some assumptions deduced 
from the observations had to be made. For instance, the 
modeled rotation curve they use and that comes from 
the bulk velocities observed along the photometric ma- 
jor axis is flat all through the inner disc (except in the 
central regions), with a maximum velocity of ^ 250 km 
s~^. Another 3D model derived from the same WSRT 
data and Hi observations of Emerson (1974) is presented 
in Braun (1991). It describes the gas and velocity struc- 
tures in terms of spiral density waves and shows a central 
disc which is tilted by ^15° from the median plane of 
the galaxy. Modeled position velocity diagrams reveal a 
complex velocity structure (sec Fig. 5 of Braun 1991). 

A complete 3D analysis of the DRAO datacube tak- 
ing into account both a warped and flaring disc, with 
spiral and/or other density waves, in addition to other 
processes like e.g. a lagging halo, as observed in other 
galaxies (Fraternali et al. 2001, 2002; Barbieri et al. 2005; 
Oosterloo et al. 2007), is beyond the scope of this arti- 
cle that aims at presenting preliminary dynamical results 
from more simple geometrical and kinematical hypothe- 
sis. We only analyze the kinematics of the bulk rotation 
of M31 by fitting a tilted-ring model to a velocity field, 
similarly to what is usually done in many other extra- 
galactic studies. The inclination of the disc is very well 
suited for such a model because the degeneracy between 
the rotation velocity and the inclination is small around 
75° (Begeman 1989). 

Probably the major ambiguity with our analysis comes 
from that it seems diSicult to decide for pixels that ex- 
hibit more than one peak which of the spectral compo- 
nents is the best tracer of the bulk disc rotation. Sev- 
eral methods can be applied to sort the different lines 
and build a useful velocity map. For instance, lines can 
be sorted by their amplitudes, integrated fluxes, velocity 
centroids or widths. There is generally one line whose 
peak amplitude or integrated intensity strongly domi- 
nates the spectrum. However choosing the brightest peak 
for the bulk rotation such as what is shown in Fig. 6 
remains problematic in pixels where the Hi emission is 
dominated by other structures than the bulk disc (like 
extraplanar gas) or when the projection effects contami- 
nates the emission. In particular, the Hi emission in the 
central regions is dominated by gas from the warped part 
of the disc. This is illustrated in the position velocity plot 
made along the major axis of the galaxy (Fig. 7; see also 
Brinks & Burton 1984). The steep innermost gradient of 
the faint central ring-like structure is dominated by the 



shallow linear emission from the warped gas. Choosing 
this brightest later line would lead to discontinuities in 
the velocity field and thus to an erroneous inner rotation 
curve. For each pixel a good compromise is to select the 
component that has the largest velocity relative to the 
systemic velocity of the galaxy while avoiding isolated 
faint features. It permits rejection of the bright compo- 
nent from that shallow linear emission, which is closer 
to the systemic velocity, and the components beyond the 
flat part of the position-velocity diagram, which come 
from possible extraplanar gas (e.g. high velocity clouds). 
The velocity field that is generated following that method 
is shown in Fig. 8 and a velocity cut along its major axis 
is reported in Fig. 7 as a green line. That map prob- 
ably allows to obtain the best velocity continuity with 
no 3D model of the datacube. Hereafter we refer to the 
"main component" the Hi line that has been chosen to 
derive this velocity field. The mass attached to the main 
component is 2.49 x 10^ Mq. 

5.1. Tilted-ring model 

The kinematical parameters and rotation curve are de- 
termined by fitting a tilted-ring model to the velocity 
field following the procedure described in e.g. Verhei- 
jen & Sancisi (2001) or Chemin et al. (2006). Briefly, 
the rotcur task (Begeman 1989) of the GIPSY data re- 
duction software (van der Hulst et al. 1992) is used to 
fit the variation of kinematical parameters as a func- 
tion of radius. It is considered here that gas motions 
are along circular orbits. Hence, no ax;isymmetric radial 
motions (gas inflow or outflow) or vertical motions per- 
pendicular to the galaxy plane are derived so that sky 
plane velocity Wobs can be written as the projection of 
the only rotation velocity along the line-of-sight Vohs = 
Vsys + Vrotcos((?) sin(i). Here 9 is the azimuthal angle in 
the galactic plane, i the inclination and Vgys the systemic 
velocity of M31. 

In the innermost regions of the disc the largest radial 
velocities are mostly not detected around the (photo- 
metric) major axis but at |F| > 0. As a consequence we 
decided to use the full velocity field to do the tilted-ring 
model analysis. We have verified that the basic results 
are not affected by the choice of the opening angle around 
the major axis by comparing with results obtained using 
a smaller angle (half-sector of 30° instead of 90°). As 
expected, the parameters derived using the full coverage 
of the kinematics are better constrained and the rotation 
curve is probed more deeply towards the galaxy centre 
for the large opening angle than for the small one. The 
shape of the rotation curve is unaffected by the opening 
angle. 

Notice however that a | cos(6')[ weight is applied to the 
data points during the fitting, giving less weight to pix- 
els close to the minor axis. This is because the contribu- 
tion of Viot to the line-of-sight velocities is less important 
along the minor axis (see also Begeman 1989; de Blok et 
al. 2008). 

One has to notice that the NE spur "N2'' and the exter- 
nal arm have been masked for the tilted-ring model anal- 
ysis because their kinematical properties make them not 
linked to the disc or the other adjacent spur-like struc- 
ture (see Figs. 5, 14 and §7). Leaving them in the velocity 
map would add large scatter in the results as well as a 
fail of the tilted-ring model at some outer radii. 
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Fig. 6. — Hi integrated emission and velocity maps of M31. The Hi components are sorted by decreasing brightness from left to right. 



The location of the dynamical centre and systemic ve- 
locity Vsys are first fitted to the velocity field by keep- 
ing the disc inclination i and position angle of the kine- 
matical major axis P.A. fixed at the photometric values 
75° and ^ 38° respectively). The fitted systemic ve- 
locity of M31 is Vsys — —304.5 i 6.8 km s ^ , whose value 
is in very good agreement with what is found using the 
integrated spectrum (§4). The location of the dynamical 
centre is found to be offset from the photometric centre 
by ~ 0.55 kpc (in projection). However, the dispersion 
around that location is ~ 0.68 kpc (in projection), show- 
ing that this offset is not really significant. 

Then the variation of i, P.A. and ?;rot are measured. 
The position angle is generally very well constrained with 
small uncertainties (see top panel of Fig. 9). It can be 
fixed in a next step to fit the inclination and the rotation 
curve. The radial profiles of i and P.A. sometimes display 
ring to ring wiggles that may look artificial (Fig. 9). This 
is the reason why the smoothed profiles are used to derive 
the rotation curve in a final step. Model velocity and 
residual velocity maps are then generated. The fitting 
of the parameters and rotation curve is repeated until 
a minimum is found for the average and scatter of the 
residual map. 

5.2. Major axis and inclination variations 

Figure 9 displays the observed and smoothed profiles 
of P.A. and i and Table 5 lists the parameters. 

We identify five distinct regions in the profiles. First 
there is a central region R < 27' where the inclination 
(position angle) strongly increases as a function of radius 

-GOO 




Angular Qffsel, 

Fig. 7. — Position velocity plot made along the photometric 
major axis P.A. = 38°. A green curve represents a slice through 
the velocity field of the main Hi component. 




-100 



300 



^-300 



3 



50" 45"^ 35" 
Right Ascension 



30" 



-400 



-500 



-600 



Fig. 8. — Velocity field of M31 used for the titled-ring model 
analysis. 

by ^ 40° (10°, respectively). Then a bump is observed 
between 27' < R ^ 62' for the position angle which 
reaches its maximum at i? = 45' {P.A.= 39°). Then an 
extended region between 62' < R ^ 120' where a very 
small increase of position angle is detected while the in- 
clination remains remarkably constant around 74°. Be- 
tween 120' < R ^ 145' a dip both in the inclination and 
position angle profiles is detected. The inclination drops 
down to 61° and the position angle to 29°. Notice that 
the uncertainties are larger in that radial range. Finally, 
beyond R = 145', the position angle remains constant 
while the inclination is observed to slightly increase to- 
wards larger values from 77° to 80°. The average Hi disc 
inclination is (74.3 ± 1.1)° and the average Hi disc posi- 
tion angle is (37.7±0.9)°, as measured for 27' < R < 120' 
(from 6 kpc to 27 kpc). In that radial range, the kine- 
matical inclination is therefore by 3° lower than the one 
derived with optical surface photometry. 

5.3. The HI rotation curve 

The Hi rotation curve of M31 is shown in Fig. 10. Its 
velocities are listed in Table 5. The rotation curve is very 
peaked in the innermost regions, showing velocities up to 
Vrot— 340 km s^^. Here, the asymmetry between both 
disc halves is prominent. A velocity dip is then observed 
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Fig. 9. — Hi disc inclination and position angle of the kinematical 
major axis as determined from tilted-ring models to the Hi velocity 
field of Messier 31. Dashed curves represent the smoothed models 
of i and P. A. used to fit the rotation velocities. Horizontal dotted 
lines show the photometric values. 
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Fig. 10. — Hi rotation curve of Messier 31. Filled diamonds 
are for both halves of the disc fitted simultaneously while blue 
downward/red upward triangles are for the approaching/receding 
sides fitted separately (respectively). 

at i? = 18.3' (4 kpc). The rotation curve is then observed 
to increase up to Wiot= 267 km s~^ at i? 47', to remain 
roughly fiat between 47' < i? < 75' (u,.ot= (264 ± 4) 
km s~^), to decrease down to ^ 230 km s"^, to remain 
flat between 95' < i? < 120' (w,.ot= (230 ± 4) km s^i) 
and finally to increase up to Wrot= 275 km s^^ in the 
outermost regions. 

The uncertainties on the rotation velocities are derived 
as follows. If gas rotation is made on purely axisymmet- 
ric and circular orbits, then rotation velocities are ex- 



actly the same at opposite sides of the galaxy. If the 
disc is perturbed, then the axisymmetry is usually bro- 
ken and differences between rotation velocities of the ap- 
proaching and receding halves exist. We thus choose a 
definition of velocity uncertainties that also takes into ac- 
count this effect. Uncertainties are indeed defined as the 
quadratic sum of the formal la uncertainty from rotcur 
with the maximum velocity difference between rotation 
curves (both disc sides model values derived simulta- 
neously minus the approaching or receding side model 
values derived separately, weighted by the number of 
points in each side). With this definition, each uncer- 
tainty is generally larger and more conservative than the 
very small one derived with rotcur. It draws a better 
representation of any departure from axisymmetry. As a 
result the largest uncertainties are observed in the cen- 
tral regions. Elsewhere the uncertainties on the rotation 
velocities are very small, generally lower than 10 km s~^. 
It shows the robustness of the method used to select lines 
to create the velocity field. Notice that for some of the 
annulii beyond R — 148' velocities could only be derived 
for the receding half. At those radii, the uncertainty is 
fixed at the formal la error provided by rotcur. 

5.4. Comparison with other works 
5.4.1. The rotation curve 

Figure 11 displays several previous Hi rotation curves 
of M31 (Newton & Emerson 1977; Brinks & Burton 1984; 
Braun 1991; Carignan et al. 2006) and our new result. 
The agreement with our previous result (Carignan et al. 
2006) is not good between 15 kpc and 23 kpc. This 
difference can be explained by the fact that Carignan 
et al. (2006) used a single emission line approximation 
while deriving the velocity field from low resolution Hi 
data of Unwin (1983). In the external parts the curves 
are in better agreement. The curve of Carignan et al. 
(2006) is flat at these radii, whose feature is not observed 
anymore. This is likely due to the assumption of constant 
position angle and inclination in our previous study. The 
agreement with the data points from Newton & Emerson 
(1977) looks better. Except for i? < 3 kpc and for some 
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Fig. 11. — Comparison between Hi rotation curves from New- 
ton & Emerson (1977, crosses). Brinks & Burton (1984, filled cir- 
cles) Braun (1991, open triangles), Carignan et al. (2006, filled 
diamonds) and our new derived rotation curve (solid line). 

points between 7 < i? < 16 kpc, the agreement with the 
curve of Braun (1991) is not good. The low velocities 
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he measures around 15 kpe and at large radius (down to 
Vrot— 190 km s^^) cannot be reproduced from the tihcd- 
ring analysis. The origin of this discrepancy may be due 
to different choices of velocity components to extract the 
velocity fields. The agreement with the model of Brinks 
& Burton (1984) is good within 4 < i? < 15 kpc and 
poor beyond 19 kpc because the rotation curve is not 
flat throughout the whole disc. 

5.4.2. The warp of MSI 

Results shown in Figure 9 imply the presence of two 
warps in M31. The first warp is located in the central 
parts (J? < 6 kpc) where the disc orientation appears 
less inclined by ~10° (in average) than the disc. The dip 
at 4 kpc thus corresponds to perturbed motions located 
in the inner warped disc. The nuclear region (R < 2 
kpc) is even more perturbed as it behaves like a "warp 
in the warp". Here rotation velocities, inclination and 
position angles are very different than anywhere else in 
the disc. Gas is observed to be very close to face-on. The 
very large uncertainties on the rotation velocities show 
that noncircular motions may be very important here. A 
second warp is detected in the outer parts {R > 27 kpc), 
where the disc is surprisingly oriented like in its inner 
warped part {P. A. and inclination) and then becomes 
more and more inclined at the largest radii. The fact 
that Hi gas is observed to start rotating faster where 
the outer warp appears is also probably a consequence of 
the inferred perturbed orbits. Notice however that the 
symmetry is excellent between both halves whereas it is 
not the case in the inner warp. 

This is not the first time that warps are evidenced in 
M31. One can mention results by CiarduUo et al. (1988) 
and Braun (1991) for the inner regions, and Newton & 
Emerson (1977), Henderson (1979) or Brinks & Burton 
(1984) for the outer regions. 

• Braun (1991) showed that the Hi distribution is 
tilted by 15° in the inner 10' (2.3 kpc). We con- 
firm the presence of the inner Hi warp in the DRAO 
observations. It is actually observed to extend up 
to i? = 6 kpc for a maximum tilting angle of ^ 15° 
and a maximum twisting angle of ^ 7° with re- 
spect to the average disc inclination and position 
angles as given above. If the perturbation in the 
nuclear region is genuine, then the warp is even 
more prominent, with maximum tilting and twist- 
ing angles of 42° and 10°. 

• The ionized gas distribution is more circular in the 
bulge, implying a more face-on disc (CiarduUo et 
al. 1988). Our measurement confirms this trend. 

• Newton & Emerson (1977) derived warp parame- 
ters in the outer disc {R > 28.5 kpc, or 125'), with 

a decreasing position angle (down to 30°) and an 
increasing inclination (up to 83°). The shape, am- 
plitude and locaion of the outer warp detected in 
the DRAO data are in agreement with their early 
result, though the inclination drops before becom- 
ing larger. 

• Another modelization was done by Henderson 
(1979), in which the warp starts at 18.2 kpc (80') 
with an increasing inclination in step of 0.75° per 



kpc, and a warped region rotated by 10° with re- 
spect to the central disc position angle. This result 
is not totally consistent with our new result. 

• In the model of Brinks & Burton (1984), the disc 
starts to warp at ~ 18 kpc, with a maximum warp 
angle of 15°. We detect a gradient of inclination 
in agreement with that value, though it is seen to 
occur at larger radius. 

According to Brinks & Burton (1984), 39% of the total 
Hi mass resides in the warped part of the disc. This 
number has to be compared with the gas mass of the 
main component integrated within R > 27 kpc (0.14 x 
10^ Mq) added to the total mass of all other spectral 
components than the main one (1.74 x 10® -Mq). It 
corresponds to 44% of the total Hi mass. This fraction 
is comparable with the prediction of Brinks & Burton 
(1984), though being slightly higher. 

6. ORIGIN(S) OF THE OTHER SPECTRAL COMPONENTS 

No kinematical analysis is presented for the other spec- 
tral components than the main disc component. Only a 
brief discussion on their possible origin(s) is proposed 
here. So many Hi peaks like those observed in M31 are 
rarely evidenced in extragalactic sources. The sum of all 
of the integrated emission other than the main compo- 
nent represents 41% of the total Hi mass. A significant 
fraction of it is due to the overlap emission from the warp 
whose signature is more prominent at large Y than along 
the major axis due to the inclination of the disc. An- 
other part is due to an additional gas component whose 
origin(s) can be manifold. 

A first origin could be structures lying in the disc it- 
self, like e.g. unresolved spiral arms by projection effects. 
Indeed, the 3'''^, 4**^ and 5^^ maps of Figure 6 show that 
projection effects may play an important role in creating 
multiple components due to the presence of more pix- 
els at \Y\ > than along the major axis. Moreover the 
projected spatial distribution of the gas in the multiple 
peaks is principally concentrated within high surface den- 
sity regions, like the Hi "ring" aroimd 13 kpc and along 
the spiral arms. Notice here that more components are 
detected in the north-western front side of the galaxy 
than on the far side. Another internal origin of multi- 
ple components could be caused by expanding gas shells 
induced by stellar winds in star forming regions. Such 
a phenomenon has already been seen in other galaxies 
(e.g. Hunter, Elmegreen & van Woerden 2001). Several 
spectra of M31 show more or less symmetric peaks cen- 
tered about the main component and could point out gas 
outflow in high density regions of the M31 disc. 

Other origins could be extraplanar gas in the form of 
e.g. a lagging halo or high velocity clouds. In recent deep 
Hi observations a lagging halo corresponds to the thick 
gaseous layer that is observed to rotate more slowly than 
the host equatorial thin Hi disc of spiral galaxies (Fra- 
ternah et al. 2001, 2002; Barbieri et al. 2005; Oosterloo 
et al. 2007). The halo emission which is often referred to 
as an "anomalous" emission has a mass that can reach 
~ 30% of the total gas mass (Oosterloo et al. 2007). It 
is thought to be gas infalling onto the disc, mainly due 
to a galactic fountain mechanism but also to the accre- 
tion from the intergalactic medium (Fraternali & Binney 
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2006, 2008). Models of gas in the halo show it has a 
larger velocity dispersion than gas in the cold disc (Oost- 
erloo et al. 2007). The DRAO datacube exhibit many Hi 
peaks whose radial velocity is observed to be closer to the 
systemic velocity than the main Hi disc component and 
whose velocity dispersion is larger than the main disc 
component. It implies that M31 could host a lagging 
halo as well. 

High or intermediate velocity clouds (HIVC) similar 
to those orbiting around the Milky Way (Wakker et al. 
1999) are observed in a few galaxies among which M31 
is a very good candidate (Thilkcr et al. 2004; Westmcicr, 
Braun '& Thilkcr 2005). M31 HIVC detected by West- 
meier, Braun & Thilkcr (2005) have a typical mass oc 10^ 
M.Q for a size of ~ 1 kpc. Similar high or intermediate 
velocity clouds could be detected in the DRAO field-of- 
view. Indeed the minimum detectable column density of 
the data corresponds to 1.6 x 10* M-q per spatial reso- 
lution element (synthesized beam size). 

One finally notices that no obvious forbidden velocities 
are detected in the high resolution datacube of M31 at 
a 3(7 detection level. Forbidden velocity clouds (FVC) 
correspond to apparent counterrotating gas, like in NGC 
2403, NGC 891 or NGC 6946 (Fraternah et al. 2002; 
Oosterloo et al. 2007; Boomsma et al. 2008). The pres- 
ence of FVC could point out ongoing accretion onto a 
host disc. As a consequence if any Hi gas accretion re- 
sponsible of part of the multiple lines is occuring in M31, 
it does not seem to be done from apparent counterrotat- 
ing material. A more careful search of HIVC or FVC in 
our datacube will require to filter the high resolution dat- 
acube in the spatial and spectral dimensions to increase 
the signal-to-noise ratio. 

7. THE Hi outskirts 

7.1. The NE and SW spurs 

Neutral gas in M31 is detected out to ^ 34.1 kpc in 
projection 36.3 kpc) in the SW approaching side (NE 
receding side respectively) . Two extended thin spurs are 
observed at the north-eastern extremity of the disc while 
gas seems more confined in its south-western extremity. 
These structures are connected to the spiral arms as seen 
in Figs. 4, 5 and 12. 

A column density of 6.2 x 10^° cm~^ is reached in the 
south-western Hi extension. A velocity gradient is de- 
tected from —515 km s~^ to —530 km s~^ for a length of 
^ 30' (6.8 kpc), in apparent continuity of the disc kine- 
matics. The velocity dispersion of the Hi peaks is 12 km 
s~^ in its northern part while and lower values are seen 
in its southern part. 

The NE spur (Nl) is oriented North-South and has a 
length of ~ 40' (9.1 kpc). Its kinematics is in continuity 
with that of the disc. A velocity gradient from —85 km 
s~^ to —65 km s~^ is detected along it while the velocity 
dispersion of the Hi lines remains constant (~ 8 — 10 
km s~^). The highest Hi density column observed in 
it is 4.8 X 10^" cm~^. A small "hook-like" structure at 
^2000 = 00'M9">11.86%(52ooo = +43'^06'24" is observed 
as part of the spur. The other spur-like structure (N2) 
is shorter (~ 30'). A velocity gradient is also observed 
along it, from —160 km s~^ to —140 km s~^. Its velocity 
dispersion is very low (~ 5 km s""'^) and is also seen to be 
constant. An average column density of 6 x 10^'^ cm^^ is 
derived for it. Notice that the three Hi extensions have 



very few pixels with multiple spectral components. 

Are these structures real? What could be their ori- 
gin(s)? One has to notice first that these extensions are 
genuine structures of the disc because they are also ob- 
served in the WSRT image (Braun et al. 2009) and they 
cannot be mistaken with residuals from a bad subrac- 
tion of MW Hi due to their kinematics. Then the kine- 
matical properties of the two NE spurs strongly differ 
from each other. We argue they are part of two different 
Hi structures, though they are located in the same area 
of the field-of-view. On one hand a close inspection of 
the datacube points out that the spur N2 seems to be a 
kinematical extension of the external spiral arm of M31 
(Fig. 5). On another hand the brightness and kinemati- 
cal properties of the spur Nl and the SW extension are 
quite similar. In particular they are responsible of the 
velocity rise that occurs at the outermost radii of the ro- 
tation curve. What is also remarkable with them is that 
they both overlay or point towards two diffuse stellar 
clumps, as seen in Fig. 13. This image displays the Hi 
column density contours superimposed onto the stellar 
distribution, as measured from deep photometry (Ibata 
et al. 2001, 2005, 2007). The SW extension is coincident 
with the Gl clump while the NE spur coincides with part 
of the northeast clump. According to Ibata et al. (2005) 
stars in the Gl clump rotates by 68 km s~^ faster than 
the disc, which is considered to rotate at ^ 230 km s^"'^ 
in their model. The true rotation velocity in the south- 
western extension is not observed constant at large radii, 
but increasing. This increase implies a lagging velocity 
of up to 45 km s~^ larger than the velocity of ~ 230 km 
. The lagging velocity of gas is thus lower than that 
of stars. The spur Nl to the North-East is observed to 
rotate faster than Wioi = 230 km s""*^ as well by a similar 
amplitude than in the southern spur, while the NE stel- 
lar clump is observed to rotate slower. Hence gas and 
stars follow the same velocity trend in the SW region of 
the disc, while it does not seem to be the case in the NE 
region. Notice however that stellar velocities are less cer- 
tain in the NE clump because of the contamination by 
MW stars. It nevertheless indicates a strong relationship 
between stars and gas at the outskirts of M31, at least 
in the SW region. 

The star formation history and metallicities of the 
northeast and Gl overdensities led Faria et al. (2007) and 
Richardson et al. (2008) to propose they could be mate- 
rial initially formed into the disc that have been stripped 
by tidal effects. The kinematics of the Hi spurs imply 
that they are bound to the disc. Though they seem to 
rotate faster than gas, it is likely that stars in the ex- 
ternal clumps are also bound to the disc of M31. A fit 
to the velocity field of a radial motion tij-ad sin(i) sin(6') 
in addition to the rotational one gives an average value 
Wrad = 6.2 km s^^ with a standard deviation of 7.4 km 
s-i for R > 150' {R > 35 kpc). 

By assuming that the spiral arms are trailing in the disc 
and by looking at the distribution of dust lanes in opti- 
cal images of M31, the front side of the galaxy is to the 
NW of the major axis and the rotation is done counter- 
clockwise. Therefore, at first order outflow motions are 
detected in the outskirts of M31 because of the derived 
positive radial velocities. At second order these radial 
motions could be mistaken by the presence of either ver- 
tical motions to the galaxy plane, which would thus be of 
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Fig. 12. — Hi emission, velocity and velocity dispersion maps zoomed in the region of the NE (top panel) and SW (bottom) structures. 



the order of 35 km for the observed inclination (80°, 
Fig. 9), or elliptical streaming in e.g. a m = 2 perturbing 
potential. However that later hypothesis is less uncertain 
here because no evident spiral structure is seen at those 
radii but rather thin "filamentary" Hi spurs. The real- 
ity is probably a combination of these hypotheses (pure 
outflow, z-motions and elliptical streaming). Vrad ~ 6 
km s""'^ is likely an upper limit of any possible outflow- 
ing motions. Such radial motions are consistent with the 
scenario proposed by Faria et al. (2007) or Richardson 
et al. (2008). Indeed, one expect outflow motions in the 
disc outskirts if tidal stripping is occuring. N-body mod- 
els would be helpful to firmly validate the detection of 
ongoing outflow in the outskirts of M31. 

7.2. The external arm 

The external arm is another new perturbed structure 
in M31. It has many properties different from other 
disc structures. It has no identifiable morphological and 
kinematical symmetric counterpart^ with respect to the 
galactic centre. Its kinematics significantly differs from 
radial velocities within the disc (Fig. 5). It docs not con- 
tain so many pixels with multiple Hi peaks. It is one of 
the faintest structures detected in the field-of-view with 
a brightness comparable to interarm regions. The total 
Hi mass attached to this structure is 1.1 x 10^ A^o- Its 
apparent length is 2.33° or 32 kpc. The kinematical dif- 
ference with the disc is caused by different orientation 
angles. Indeed its position angle and inclination are ~ 

® The only noteworthy Hi feature which could be a counter- 
part on the other edge of the disc is located around 02000 = 
00'^39™49.80=,<52000 = +41'*13'39". However it is not an extended 
structure like the external arm. 



10° larger than those of the disc (Fig. 14), with the ex- 
ception of the region B > 34 kpc where both inclinations 
become very similar. The rotation velocities of the ex- 
ternal arm are comparable with the disc rotation curve. 
No obvious stellar counterpart to the gaseous arm can 
be identified in the optical images of Ibata et al. (2005, 
2007) because of their too high contrast. 

Here again one may question about its origin. Inter- 
nally driven perturbations are hardly possible because 
one would expect the creation of axisymmetric or bisym- 
metric features in the gas distribution (e.g. a bar, a ring 
or two spiral arms) . The apparent isolated nature of the 
arm on one half of the disc and not on the other side rules 
out this hypothesis and points out to an external pertur- 
bation for its origin. It includes spiral arm triggering 
by tidal effects exerted by a companion or even gas ac- 
creted from a companion. This would be more consistent 
given the perturbed nature of the stellar distribution at 
large radius caused by tidal interactions with low mass 
satellites. Indeed the stellar halo of M31 is filled of stel- 
lar residuals from past interactions with smaller galaxies. 
M31 has currently two nearby bright satellites, M32 and 
NGC 205 at projected distances of ~ 25 kpc and ~ 37 
kpc to the M31 nucleus. Both are expected to have inter- 
acted with M31. In a gas accretion scenario, one would 
need an efficient tidal stripping because both companions 
are devoid of gas. Block et al. (2006) simulated a head-on 
collision with M32 that could have occured 210 million 
years ago and have generated the ring-like morphology 
in the distribution of dust and gas. Howley et al. (2008) 
simulated an interaction with NGC 205 and found radial 
orbits at high velocity (up to 500 km s~^) for the small 
companion. According to the possible orbits from the 
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Fig. 13. — Comparison of the Hi distribution (white eontours) 
with the faint extended stellar disc of M31 from Ibata et al. (2007, 
grey-scale image and grey contours adapted from their Fig. 50.). 
The letter A refers to the North-East stellar clump, B to the Gl 
clump and C to the Giant Stream, following the nomenclature of 
Ibata et al. (2007). Hi contours are for column densities of 5 X 10^^, 
1 X 10^0, 4 X 10^", 9 X 10^0 and 2 X lO^l cm'^. 
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Fig. 14. — Hi inclination and position angle of the external arm 
(open symbols). Filled symbols show the orientation of the disc 
parameters from Fig. 9. 

models, we think the most hkely impactor that gener- 
ated the external arm could be NGC 205. The location 
of the perturbed external arm indeed coincides with the 
path of NGC 205 in the apparent far side of the galaxy 
(see Fig. 19 of Howley et al. 2008) while the M32 trajec- 
tory leads the compact elliptical close to the M31 centre 
along a polar orbit. Here again more detailed numerical 



simulations are needed to investigate the different possi- 
bilities and find the origin of this perturbation. 

8. MASS DISTRIBUTION ANALYSIS 

The mass distribution modeling of the galaxy is done 
by decomposing the total gravitational potential into a 
supermassive black hole, a luminous baryonic component 
and a dark matter component. It consists in fitting a 
dynamical model for the dark matter component to the 
observed rotation curve, taking into account the baryonic 
rotation curves. Two or three free parameters are fitted 
for the halo in addition to the stellar mass-to-light ratios 
for the bulge and galactic disc. These latter can be fixed 
under some assumptions described below. 

8.1. The central supermassive black hole 

Messier 31 is known to have a central supermassive 
black hole which mass can now be robustly constrained. 
The black hole contribution to the rotation curve remains 
negligible through the disc because of the point mass 
nature of its gravitational potential. A mass of 1.4 x 
10* Mq was derived using Hubble Space Telescope STIS 
spectroscopy of the nucleus (Bender et al. 2005) and is 
used in our study. 

8.2. The luminous baryonic matter 

The contribution of the luminous baryons to the rota- 
tion curve consists in stellar and gaseous components. 

8.2.1. The stellar potential 

The gravitational potential of the stars is decomposed 
into contributions from a bulge and a disc. These ones 
are derived from the i?— band surface brightness profile 
of M31 (Walterbos & Kennicutt 1987, 1988). 



o 
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Fig. 15. — i?— band surface brightness profile of Messier 31. Ob- 
servations are from Walterbos & Kennicutt (1987). Filled symbols 
are the observations, a red dashed line is the adopted contribu- 
tion of the disc to the surface brightness, a blue dotted line the 
adopted contribution of the bulge and a solid line the sum of both 
components. A green dashed line shows the fit of the exponential 
disc. 

A bulge-disc decomposition of the surface brightness 
profile is done using an exponential disc contribution for 
the stellar disc plus a de Vaucouleurs i?^/^ law (de Vau- 
couleurs 1948) for the bulge contribution. However, the 
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do Vaucouleurs law gives a bad result for the bulge con- 
tribution to the surface brightness profile because it sig- 
nificantly overestimates the intensity in the inner parts 
of the galaxy. A more robust bulge-disc decomposition is 
done using a generalized Sersic R^^"' model for the bulge 
light (Sersic 1963, 1968) plus an exponential disc. In the 
Sersic model, the intensity profile is written as 



I{R) = le exp {-6„ [{R/ReY/'' - l] } 



(1) 



in which le is the intensity at the effective radius Re de- 
fined as the radius that contains half of the total light 
and n a dimensionless index that measures the "curva- 
ture" of the luminosity profile. 6„ is a function of n and 
is defined as 6„ = 1.9992n- 0.3271 for 0.5 < n < 10 (Ca- 
paccioli 1989). The de Vaucouleurs law is thus a specific 
case of the Sersic law for n = 4. 

The stellar disc intensity profile is fitted with the ex- 
ponential disc formula (see e.g. Freeman 1970) 

7(i?) = /o exp(-i?/i?d) (2) 

where Iq is the intensity at iZ = and Rd the scale-length 
of the disc. 

The two components were fitted simultaneously to the 
brightness profile (Fig. 15). The fitted bulge parameters 
are Rg = (1.3 ±0.1) kpc, /Xe = (18.6 ±0.1) mag arcsec"^ 
(not corrected for foreground and internal dust attenua- 
tion) and r?, = 1.1 ± 0.1. The total apparent magnitude 
of a bulge for the Sersic law is given by : 



m = He — 5 log Re — 2.5 log 



27m 



2n 



r(2n) 



(3) 



(Graham & Driver 2005), which gives a total uncorrected 
apparent magnitude of meuige ~ 3.25 mag for M31. 

The fitted disc parameters are Rj, ~ (5.6 ± 0.1) kpc 
and a central surface brightness /xq = (19.0 ± 0.1) mag 
arcsec"^ (not corrected for foreground and internal dust 
attenuation nor inclination effects). This disc scale- 
length compares well with values fitted with the same 
dataset (Walterbos & Kennicutt 1987, 1988; Geehan et 
al. 2006). As seen in Fig. 15, the adopted disc brightness 
profile is nevertheless not the fitted one, but the result of 
the subtraction of the fitted bulge contribution to the ob- 
served profile. This allows us to keep the stellar surface 
brightness as close as possible to the observed profile. 
Moreover, it enables us to keep the intensity variations 
of the disc that could reproduce wiggles in the rotation 
curve. The surface brightness profiles are then corrected 
from Galactic extinction (following Schlegel, Finkbeiner 
& Davis 1998) and for inclination effects (only for the 
disc profile) and for internal extinction as described in 
§8.2.2. 

8.2.2. The mass-to-light ratios of the stellar component 

Two free parameters of the mass models are the mass- 
to-light ratios of the bulge Tsuige and of the galactic 
disc Tdisc, which are considered constant as a function 
of radius for the remainder. Models presented here also 
enable these parameters to be fixed. A relative approxi- 
mation for them can be deduced from stellar population 
synthesis (SSP) models (e.g. Bell et al. 2003), provided 
that a color is given. 

The bulge color is B — R = 1.75 mag (Walterbos & 
Kennicutt 1988), which gives a color of B — R ^ 1.27 af- 
ter correction for foreground and internal dust extinction 



effects. Here wc have applied the reddening law for M31 
derived by Barmby et al. (2000), using a mean redden- 
ing factor Eb-v = 0.24 mag (Wilhams & Hodge 2001). 
This reddening law is relatively comparable with the one 
measured for the Milky Way. The corrected bulge mass- 
to-light ratio is Tsuigc ~ 2.2 Mq/Lq, following the pre- 
scriptions of Bell et al. (2003). These authors indicate a 
r.m.s. scatter of the order of 25% on the mass-to-light 
ratios. An independent constraint of the bulge mass- to- 
light ratio can be obtained directly from measurements 
of the bulge effective stellar velocity dispersion, with the 
assumption that the dark halo has a negligible contribu- 
tion to the velocity dispersion. The estimated dynamical 
mass of the bulge is A^Buige = 1-9 x 10^'^ Mq (Marconi 
& Hunt 2003), which value becomes A^suigo = 2.4 x 10^° 
Mq after scaling to our adopted bulge effective radius. 
It implies a corrected dynamical mass-to-light ratio ratio 
of Tsuige ~ 0.8 Mq/Lq. This is about 2.8 times lower 
than the ratio derived from SSP models. 

The average disc color is i? — i? ~ 1.6 mag, as derived 
from Walterbos & Kennicutt (1987), giving a corrected 
color B ^ R ^ 1.1 mag and a corrected disc mass-to- light 
ratio Tdisc ~ 1-7 Mq/Lq. When no correction from 
internal extinction is applied to the B — R colors, SSP 
models imply Tsuige ~ 4.0 and Toisc ~ 3.2. Our val- 
ues of corrected and uncorrected Tdisc are thus in good 
agreement with those estimated by Widrow & Dubinski 
(2005) or Geehan et al. (2006). 

Applying the same foreground and internal corrections 
as above (using Ab = 0.67 mag, de Vaucouleurs et al. 
1991), the total bulge luminosity of M31 is Leuigo ~ 
2.9 X 10^" Lq in the i?— band, which translates into a 
total corrected bulge mass of A^Buige ~ 6.4 x 10^° A4q 
for the adopted Tsuige = 2.2. The deduced total lu- 
minosity of the disc (also corrected for inclination) is 
Loisc ^ 4.2 X 10^" Lq, corresponding to a total cor- 
rected stellar disc mass of AIdisc ^ 7.1 x 10^" A4q for 
the adopted Tdisc = 1-7. The total inferred stellar mass 
is Mi, 1.35 x 10^^ Mq. SSP models thus infer disc 
and bulge masses comparable with each other, which 
is somewhat unexpected for a high surface brightness 
galaxy (Courteau & Rix 1999). If one uses the dynami- 
cal value Tsuige = 0.8 instead of the SSP value, then the 
total stellar mass becomes A^* ~ 9.5 x 10^" Al©, making 
the bulge-to-disc mass ratio falling down to 34%, which 
seems more realistic. 

8.2.3. Infrared surface photometry 

Determining mass-to-light ratios for galaxies in the 3.6 
/Ltm band from stellar synthesis models is not an easy task 
because no straightforward recipe exists at this wave- 
length, contrary to optical bands. 

We have attempted to derive stellar mass-to-light ra- 
tios and masses from observations performed with the 
SPITZER-lnb&ved Array Camera at 3.6 fim (Barmby et 
al. 2006). However no realistic stellar masses and mass- 
to-light ratios could be estimated, using either the for- 
malism developed in de Blok et al. (2008) that derives 
mass-to-light ratios from the Large Galaxy Atlas of the 
Two Micron All Sky Survey (2MASS) J-K bulge and 
disc colors, or the average values of bulge and disc mass- 
to-light ratios inferred for seven barred and unbarred 
Sab, Sb and Sbc spirals from the THINGS galaxy sample 
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Fig. 16. Hi surface density profile of Messier 31. A full line 
is for the total Hi emission, a dashed line for the emission of the 
main component. 



of de Blok et al. (2008). Such values both lead to either 
a too large total stellar mass or a too massive bulge com- 
pared with the disc. As a consequence mass distribution 
models have been fitted using the only ii— band photom- 
etry for clarity and simplicity reasons. 

8.2.4. The gaseous disc potential 

The contribution to the rotation curve from the 
gaseous component is derived from the H2 and Hi sur- 
face densities. The H2 surface density is taken from the 
CO gas survey of Dame et al. (1993), after scaling to our 
adopted distance. The total molecular mass is ~ 3.5 x 10^ 
M.Q. The Hi surface density profile is determined by av- 
eraging the total Hi emission over elliptical rings which 
orientations are given by the tilted-ring model described 
in section 5.1. 

The Hi surface density profile is displayed in Figure 16, 
for the total Hi distribution and for the principal emis- 
sion component (the one that served to build the velocity 
field). The profile of the total Hi distribution is listed in 
Tab. 5. It has been derived using orientation parameters 
of Fig. 9 and Tab. 5. No correction due to the overlap 
from the warp emission has been applied to the density 
profile because no 3D model of the datacube has been 
attempted in this article. It has no consequences on the 
following result because of the very low induced ampli- 
tudes of the atomic gas rotation velocities. 

The majority of Hi gas is located between R = 7 kpc 
and i? = 18 kpc, which corresponds to the location of the 
brightest "rings". Peaks are also observed in the inner 
region, in agreement with the locations of other inner 
ring-like structures. The rotation curve of the atomic 
gas component is then derived after scaling the density 
profile by a factor of ^ 1.3 in order to take into account 
the helium contribution. The contribution of the atomic 
gas to the overall rotation curve is very low and does not 
exceed ~ 70 km s~^. As the mass of the molecular gas is 
only ^ 8% that of the total atomic gas, its contribution 
to the rotation curve is also very low (~ 15 km s^^ at 
maximum) and comparable with the velocities due to the 
black hole at the largest distances. The mass-to-light 
ratios of the molecular and atomic components are kept 
fixed during the fittings (Tgas = 1)- 



8.3. The dark matter halo potential 

Three different halos are fitted in this study. A first 

model is associated with a pseudo-isothermal sphere 
while two other models are cosmological halos as their 
density profile is derived from Cold Dark Matter (CDM) 
numerical simulations. No attempts to fit the mass dis- 
tribution in triaxial models were done as all models con- 
sider a spherical halo. Also, no attempt to model an 
adiabatic contraction of the dark halo in its central parts 
was done (Dutton et al. 2005). 

8.3.1. Navarro-Frenk-White halo 

Navarro, Prenk & White (1996, 1997) were among the 
first ones to propose a formalism that was fitted to re- 
sults from numerical simulations done in the framework 
of the CDM theory. They found a halo shape that is 
independent of the halo mass, which is why this halo is 
often referred to as the "universal" halo. The mass den- 
sity profile is steep (cuspy) as it scales with at low 
radius and is written as 



Pnfw(-R) = 



Pcrit 



{R/Rs){l + {R/Rsy 



(4) 



where pcrit = 3i7Q / (87rG') is the critical density for clo- 
sure of the Universe, Sc represents a characteristic den- 
sity contrast and Rg a scale radius. The circular velocity 
profile corresponding to this halo allows to fit two pa- 
rameters, V200, the velocity at a the virial radius -R2005 
and c — R2ml Rs, a concentration parameter of the halo, 
to the observed rotation curve, and is written as : 



V{R) = 



2001 



1 ln(l + cx) — cx/{l + cx) 
X \n{\ + c)-c/{l + c) 



(5) 



In this equation, x = R/R2oo- This model is referred to 
as the "NFW" model or NFW halo hereafter. 

8.3.2. Etnasto halo 

Merritt et al. (2005, 2006) used more recent ACDM 
numerical simulations to model the mass density profile 
of dark halos and found that they can be fitted with a 
three-parameter model. The density is written as 



Pci„(i?) = PcCXp {-d„ [(i?/i?e)'/" - 1] } 



(6) 



This relation is very similar to the expression of the lu- 
minosity profile for elliptical galaxies and bulges of galax- 
ies (Eq. 1), with the difference that the Sersic formula ap- 
plies to the (projected) surface distribution of light from 
galaxies whereas Eq. G applies to the spatial distribution 
of the halo mass. Like in Eq. 1, n is a dimensionless pa- 
rameter which measures the "curvature" of the density 
profile, dn is a function of n and can be approximated by 
dn ~ 3n- l/3-H0.0079/n (Merritt et al. 2006), provided 
that n > 0.5. Here, pc is the density at the effective ra- 
dius Re defined as the radius of a sphere in which half of 
the total halo mass is contained. 

The mass profile of the Einasto halo (Cardone, Piedi- 
palumbo & Tortora 2005; Mamon & Lokas 2005; Merritt 
et al. 2006) is written as 



Mei„(i?) = 47rnii^pee''"rf„-^"7(3n,a;). 



(7) 
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Here 7(3n, x) is the incomplete gamma function de- 
fined as ^ 

■y{3n,x)= [ e-V"-Mt. (8) 
Jo 

Merritt et al. (2006) and Graham et al. (2006) refer 
to this model as the Einasto halo due to the early works 
made by Einasto (1965, 1968, 1969) and Einasto & Haud 
(1989) on the light and mass distributions of galaxies, 
independently from those of Scrsic (1968). Simulated 
empirical galaxy-sized halos have typical values n = 5 — 7 
(Tab.l of Merritt et al. 2006), log{pe) ~ -5.5 {Mq pc'^) 
and Te of several hundreds of kpc. Hereafter, this model 
is referred to as model "EIN" or Einasto halo. 

8.3.3. Core-dominated halo 

A halo having a volumic mass density which is constant 
at low radii is often referred to as a core dominated halo 
or a (pscudo-)isothcrmal sphere. Several formalisms exist 
to describe such a halo shape (Blais-Ouellette, Amram & 
Carignan 2001, and references therein) and the following 
prescription is used : 

Piso{R) = T^, (9) 

where po and Rc arc the central mass density and the 
core radius of the halo. These are the parameters that 
are fitted to the rotation curve. This dark halo was used 
in Carignan ct al. (2006). The volumic mass density 
decreases at large radius like p oc , exactly as for the 
NEW halo. We refer to this halo as the "ISO model" 
and pseudo-isothermal or core halo. 

8.4. Fittings 

Levenberg-Marquardt least-squares fittings to the ro- 
tation curve of M31 are done. One model uses Tdisc and 
f Bulge fixed at the SSP values of 1.7 and 2.2 (model re- 
ferred to as "SSP"). Another hybrid model ("HYB") is 
fitted with Tdisc = 1-7 and Tsuigc = 0.8 (sec §8.2.2). 

We have also considered best- fit models ("BE") for 
both photometric bands, which have free halo and stellar 
parameters. Only one plot is shown for this model given 
that they lead to either very unconstrained parameters 
(e.g. the core radius), or too massive stellar disc (the 
so-called "maximum disc" fitting) or unphysical models 
(e.g. null mass-to-light ratios). This later result is not 
acceptable for a high surface brightness galaxy fike M31 
and has to be rejected. Results of best-fit models are 
listed in Tab. 4 and are only used to provide a very up- 
per limit of the total stellar mass of M31. 

A usual normal weight is applied to the data points 
as the inverse of the square of the velocity uncertainties. 
Several tries are done with different sets of initial guesses 
in order to avoid as much as possible local minima. Fur- 
thermore, a constraint is applied to the mass densities 
which have to be positive in order to avoid hollow ha- 
los. Another final constraint is for the quantity n of the 
Einasto halo, which has to be larger than or equal to 0.5 
because of the definition of d„. 

8.5. Results 

Rotation curve decompositions are displayed in Eig- 
ure 17 and fitted parameters listed in Tab. 4. 
The basic results can be summarized as: 



• The reduced are high, it is impossible to re- 
produce the exact perturbed shape of the rotation 
curve. 

• The rotation velocities in the nucleus are not repro- 
duced. It is expected since this region has a very 
negligible weight in all fittings because of the large 
observed uncertainties. 

• All models with a fixed Teuigc = 2.2 significantly 
overestimate the dip at 4 kpc, sometimes with a 
velocity difi'erence of up to 60 km s~^, as seen in 
the top-left panel of Fig. 17. 

• Hybrid models with a fixed Teuige = 0.8 gives bet- 
ter results than pure SSP models. 

• For each of the HYB and SSP model, the quality 
of the fit is equivalent whatever the shape of the 
dark matter halos is. 

8.6. Analysis 

High Xr values are in majority caused by the pecu- 
liar shape of the high resolution rotation curve, which is 
hard to reproduce irrespective of the halos and mass-to- 
light ratios used during the fittings. From a statistical 
point-of-view, the derived Xr combined with about 95 de- 
grees of freedom indicate that the difference between the 
rotation curve and each dynamical model is extremely 
significant, so that it could be tempting to reject any 
mass distribution models. However the purpose of this 
discussion is not to claim that the ultimate dark matter 
halo of M31 has been found, but rather to put basic con- 
straints on its parameters and to provide a summary of 
the dynamical content of the galaxy. 

8.6.1. Mass-to-light ratios 

We first discuss the results from SSP and HYB mod- 
els. The worst results from all purely SSP models are 
evidence that a color-based mass-to-light of 2.2 for the 
bulge has to be ruled out by the present observations un- 
less the measured rotation curve is completely uncorrect 
within i? = 5 kpc , which is hard to imagine. The bulge 
is indeed too massive in the SSP model, as seen in the 
top- left graph of Fig. 17. 

The bulge in hybrid models still overestimates the 4 
kpc dip, but the difference is not as significant as for 
the SSP models. Are stellar synthesis models compatible 
with Tsuige = 0.8? Reconciling them with Tsuigc — 0.8 
would require a dust-free color B — R^ 0.62 mag and an 
intrinsic reddening Eb-r ^ 1 mag towards the direction 
of the M31 bulge. However this value is 0.65 mag larger 
than the one adopted here. The global reddening law 
for M31 has been obtained by averaging data from many 
linc-of-sights and having different dust content (see eg 
Fig. 8 of Barmby et al. 2000). In view of the upper lim- 
its allowed by the uncertainties on Eb-r and Eb-v (as 
quoted in Barmby et al. 2000 and Williams & Hodge 
2001) and of the spatial distribution of the (largest) ob- 
served reddening factors [Eb-v > 0.4 mag), which are 
not really concentrated in the bulge of M31 but rather in 
the disc dust lanes, one can hardly argue that large red- 
denings {Eb-r > 0.7 mag) are likely in the direction of 
the bulge of M31. The adopted extinction law and galaxy 
color thus seem realistic. Even by taking into account a 
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Fig. 17. — Mass distribution models of Messier 31 
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black dashed-dotted line is for the black hole contribution, solid and dashed blue lines for the neutral and molecular gaseous discs, solid 
and dashed red lines for the stellar disc and bulge, a green line for the dark matter halo and a solid black line for the overall model. 



lower value Teuigc ^ 1-6, as allowed by the uncertainties 
on B — R and on the population synthesis models, the fit 
would remain of bad quality {Xr — 22.4 with the pseudo- 
isothermal halo for instance). Notice that Tsuige ~ 1-6 
still remains at odds with the dynamical value inferred 
from velocity dispersion measurements. We are thus left 
at first sight to question about the validity of the stellar 
population synthesis models of Bell et al. (2003) for the 
bulge of M31. We do not know the implications of a fine 
tuning of ingredients of stellar synthesis models on the 
stellar evolution of a bulge with Tsuige ~ 0.8. A further 
analysis should be conducted to study this problem. 

Adopting TBuigc = 0.8 as the most likely value implies 
a bulge mass-to-light ratio ^2.1 lower than the disc one. 
The stellar disc of M31 is massive but sub- maximum. At 
a radius R = 2.15i?d ~ 12 kpc, which is the location 
of the peak velocity of an exponential disc, the ratio of 
the stellar disc velocity to the rotation curve is ~ 45% for 
Tfoisc = 1-7. Notice that the ratio of the total bulge-|-disc 
velocity to the rotation curve rises to ^ 54%, which still 
shows that the stellar component is sub-maximum. Fig- 
ure 17 shows that the dark matter component dominates 
everywhere the mass of the stellar disc, and more impor- 
tantly for ACDM halos. 

One would need Tdisc ~ 4.8 for the galaxy to be ex- 
actly 75% of the maximum disc. This is the reason why 
the best-fit models correspond to the maximum disc hy- 
pothesis because they provide the largest mass-to-light 
ratios. There is no real concensus on the nature of discs 



to be maximum or sub-maximum since all trends are 
found in the observations and/or the numerical simula- 
tions (e.g. Bottema 1997; Courteau & Rix 1999; Weiner, 
Sellwood & Williams 2001; Kassin, de Jong & Weiner 
2006; Zanmar Sanchez et al. 2008). A point to notice 
here is that the population synthesis result seems not to 
correspond to the maximum disc solution, even though in 
general both assumptions are thought to be tightly linked 
(Bell & de Jong 2001). The upper limit of the total stel- 
lar mass implied by the BF models is 2.75 x 10^^ A4q. 
Those masses strongly differ from other stellar mass esti- 
mates of M31. Notice that BF models give either uncon- 
strained core radii or no stellar components, as seen in 
Tabs. 4. The concentration and circular velocity of the 
NFW halo are also not consistent with expectations of 
similar M31-sized halos from CDM simulations (Bullock 
et al. 2001; Neto et al. 2007). For all these reasons, the 
results of BF models are rejected in the following of the 
analysis. 

A total stellar mass of M-, = (9.5 ± 1.7) x 10^*^ Mq is 
derived using Toisc 1-7 and Tsuigo = 0.8 and a con- 
servative 25% uncertainty in the disc mass-to-light ratio. 
The comparison with results of mass distribution mod- 
elings of M31 from other works (Kent 1989; Carignan et 
al. 2006; Widrow, Perrett & Suyu 2003; Geehan et al. 
2006; Tempel, Tamm & Tenjes 2007) is limited because 
uncorrected colors or other different photometric bands 
have been used for deriving the stellar mass-to-light ra- 
tios as well as different rotation curves. An interesting 
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result raised by Widrow, Perrctt & Suyu (2003) is that 
most of their numerical models point out to a higher disc 
mass-to-light ratio than the bulge one, which is consis- 
tent with our current analysis. However Geehan et al. 
(2006) deduce a higher bulge mass-to-light ratio than for 
the disc (3.9 and 3.3 respectively), as well as Tempel, 
Tamm & Tcnjcs (2007). The total stellar mass compares 
very well with those inferred by Widrow, Perrett & Suyu 
(2003), Geehan et al. (2006) or Barmby et al. (2007), 
but is higher than in Tempel, Tamm & Tenjes (2007, 
~ 7.5 X 10^° Mq). Widrow, Perrett & Suyu (2003) ar- 
rive at the conclusion that the disc mass is < 8 x 10^" 
Mq and the bulge mass is 2.4 x 10^° A^©, Geehan et al. 
(2006) report 7.2 x 10^° Mq for the disc and 3.2 x 10^° 
Mq for the bulge. The disc mass we deduce agrees with 
their results within the uncertainties while the bulge only 
agrees well with Widrow, Perrctt & Suyu (2003). It is 
worth notifying that the stellar mass quoted in Carignan 
et al. (2006, 2.3 x lO^^ Mq) is too high compared with 
all other estimates. This result is likely attributed to the 
too large contribution of the bulge and to the use of the 
best-fit/maximum disc hypothesis in our previous mass 
distribution model. 

8.6.2. The dark halo of M31 

The only valid hybrid models are used to discuss the 
basic parameters of the dark halo. We do not aim at 
determining an exact density shape of the dark halo of 
the Andromeda galaxy (core versus cusp). We only no- 
tice that CDM models describe slightly better the mass 
distribution than the ISO model. However this result is 
not really significant due to the high inferred x^. The 
parameters of the Einasto halo are in good agreement 
with those of empirical galaxy-sized halos modelized in 
numerical simulations (Merritt et al. 2006). The effective 
radius Rg is large and not well constrained. Actually, the 
Einasto model does not help to improve the quality of the 
fit, though it has one more free parameter than ISO or 
NEW models. 

The concentration parameter of the NEW model is c = 
20.1 ±2.0. This value is consistent with the one expected 
for M31-sized halos in numerical simulations (8 < c < 20 
Bullock et al. 2001; Ncto et al. 2007). It is 2.5 times lower 
than the one derived by Seigar, Barth & Bullock (2008), 
who used a rotation curve for the dynamical analysis as 
well. 

8.6.3. Dynamical content summary 

The dynamical content of M31 is calculated by averag- 
ing results from the hybrid model of the NEW halo. The 
total dynamical mass enclosed in the inner 38.1 kpc of 
M31 at the last measured data point of the Hi rotation 
curve is A^Dyn = (4.7±0.5) x 10^^ Mq. It represents the 
sum of the dark component and baryonic masses. It im- 
plies a dark-to-baryonic mass ratio of A^Dark/A^Baryon ~ 
3.7 (79% of dark matter), and a total dynamical mass-to- 
light ratio Aloyn/Ltot ~ 6.7. Here A^Baryon represents 
the sum of the black hole, gaseous and stellar masses, 
A^Dark thc dark matter mass and Ltot the total stellar 
luminosity in the i?— band, all values integrated within 
R = 38.1 kpc. The atomic and molecular gas to stellar 
mass fraction is A^Hi+He+H2/Al* = (6.2ti;o)%- At the 
virial radius of the NEW halo (i? = 159 kpc), the total 
enclosed mass of M31 is Alvir ~ 1-0 x 10^^ Mq. 



Many works reported estimates of the total mass of 
M31 using different observational techniques to probe 
its potential. Eor instance a value of 2.8 x 10^^ Mq 
inside 31 kpc has been derived by Evans et al. (2000) 
from planetary nebulae radial velocities. Recent single- 
dish measurements of the rotation curve out to ~ 35 
kpc led Carignan et al. (2006) to the conclusion that a 
total mass of 3.4 x 10^^ Mq was enclosed inside this ra- 
dius. These values are somewhat smaller than the new 
mass derived inside 38.1 kpc. As for its total mass Evans 
et al. (2000) have used radial velocities of fifteen dwarf 
satellites of M31 to estimate a value of l.Qt^^^^ x 10^^ 
Mq. Evans & Wilkinson (2000) indicate a mass of 
12.3^6^ X 10^^ Mq, from planetary nebulae, globular clus- 
ters and galaxy satellites radial velocities measurements 
inside a very large distance to the M31 nucleus (550 kpc). 
More recent studies of the kinematics of several substruc- 
tures located in the giant stellar halo surrounding M31 
help getting more information about the total mass of the 
galaxy. Among others, Ibata et al. (2004) measure the 
kinematics of stars inside the Giant Stream and model 
their orbits within a universal NEW profile. They de- 
duce a total mass of 7.5^^ ^ x 10^^ inside 125 kpc, for 
a lower limit of 5.4 x 10^^ Mq. With more kinematical 
data of RGB stars from the stellar halo, the same team 
indicates a higher virial mass > 9 x 10^^ Mq (Ghapman 
et al. 2006). In the same time, numerical simulations of 
the merger with the Giant Stream progenitor indicate a 
best-fit mass of (7.4 ± 1.2) x lO" Mq within 125 kpc 
(Eardal et al. 2006), very close to the result of Ibata et 
al. (2004). Analyses of compilations of a lot of kine- 
matical data (rotation curves and/or velocity dispersion 
profiles of the M31 disc) allowed Geehan et al. (2006) 
to constrain a mass of 7.1 x 10^^ Mq, 10.7 x 10^^ Mq 
inside 400 kpc for Tempel, Tamm & Tenjes (2007). Ei- 
nally, Seigar, Barth & Bullock (2008) find a virial mass 
of (7.3 ±0.2) X 10" Mq for their NEW model of another 
composite rotation curve of M31. Our new extrapolated 
virial mass is thus in good agreement with the upper lim- 
its provided by Ibata et al. (2004) and Ghapman et al. 
(2006). 

9. CONCLUSION 

A new deep Hi mapping of Messier 31 obtained at 
DRAO has been presented. It combines single dish an- 
tenna and high resolution synthesis telescope observa- 
tions. The high spectral resolution of the observations 
show evidence that the most important contamination 
by the Milky Way Hi is done by a spectral component 
at ^ —40 km s^^. The main results from the analysis 
of the neutral gas distribution, kinematics and dynamics 
are : 

• Up to five Hi peaks are detected in spectra of the 

M31 disc, though the detection of four or five peaks 
remains very rare. The majority of M31 spectra 
exhibit a maximum of two Hi lines. The origin of all 
these additional lines may be internal (projection 
effects of gas in the warped part of the disc and in 
inner unresolved spiral arms, expanding clouds in 
star forming regions) and/or external (extraplanar 
gas). 

• Ring-like and spiral-like structures are observed in 
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the gas distribution. This is consistent with previ- 
ous Hi images of M31. 

• New Hi structures arc revealed in the outskirts of 
the disc. First, thin extensions (called the Hi spurs) 
are seen to the North-East and South- West. Their 
size is 7 — 10 kpc. Evidence for velocity gradients 
in those structures are presented. Then a faint, ex- 
tended (32 kpc) and outer Hi structure (called the 
external arm) is discovered. Contrary to every Hi 
structures identified within the field-of-view, it has 
no morphological nor kinematical counterpart in 
the opposite side of the galaxy. Its peculiar kine- 
matics can be explained by a 10° larger inclina- 
tion and position angle than those of the disc. One 
of the Hi spur appears as a genuine kinematical 
extension of that external arm. 

• A tight relationship is evidenced between the Hi 
spurs and perturbed stellar structures (the "Gl" 
and north-eastern clumps) . As possible radial out- 
flowing motions are detected in the Hi spurs, we 
argue that the perturbed stellar and gaseous out- 
skirts are being torn from the M31 disc. Further- 
more the external arm could be the vestige of a 
past tidal interaction, or even perhaps gas accre- 
tion onto the M31 disc. 

• The rotation curve is peculiar. In particular it 
shows a central velocity dip at 4 kpc and a ris- 
ing shape in its outer part. A peak up to 340 km 
s~^ is observed within the inner 2 kpc. Except for 
the inner regions, the axisymmetry of the gas rota- 
tion between the two disc halves is very good. The 
rotation velocity uncertainties are mostly below 10 
km s~^. 

• Between i? = 6 kpc and R = 27 kpc a mean incli- 
nation of (74.3 ± 1.1)° is derived for the Hi disc of 
M31, for a mean position angle of (37.7 ± 0.9)°. A 
prominent Hi warp is detected inside R = 6 kpc, 
where the disc appears less inclined, while another 
external warp detected beyond R = 27 kpc makes 
the disc being less inclined and then more inclined. 
The disc tilting is accompanied by a twist of the 
kinematical major axis. 

• The peculiar shape of the rotation curve is hardly 
reproduced by mass distribution models. The cen- 
tral velocity dip cannot be modelled by any fittings. 
The bulge mass-to-light ratio deduced from color- 
based stellar synthesis models appropriate. A value 
close to the expected one from stellar velocity dis- 
persion measurements provides better results. 

• A total stellar mass of Mi, = (9.5 ± 1.7) x 10^° 
A^0 is found. A dynamical enclosed mass of 
Mi)yn = (4.7 ± 0.5) X 10" Mq is derived at the 
last measured radius of the rotation curve (i? = 38 
kpc). A dark-to-baryonic mass ratio of ~ 4 and a 
total gas-to-stellar mass fraction of (6.2li[ ;q)% are 



also derived inside that radius. The total mass of 
M31 extrapolated to its virial radius {R = 159 kpc) 
is Mvir ~ 1.0 X 10^^ Mq. Both these values are 
consistent with those given by other works using 
other dynamical tracers or methods. 
Several future works should be conducted with these 
data. A first crucial point will be to investigate the origin 
and evolution of the perturbed Hi outskirts - the external 
arm and spurs. With the help of future N-body and 
hydrodynamical simulations, the Hi kinematics should 
put new constraints on the mass assembly of M31. A 
second challenge will be to understand the origin of the 
peculiar velocities in the central 4 kpc. The perturbed 
velocities are related to the presence of the inner warp, 
which has perhaps been generated by an interaction with 
the small companion M32 (F. Combes fc F. Bournaud, 
private communication) and/or by the bar/boxy bulge 
perturbation. Another work will focus on a 3D modeling 
of the datacube in order to study the origin of the gas 
components which are seen in addition to the warped 
disc of M31. Whether M31 contains the "anomalous" gas 
from a lagging halo, as it is observed in recent deep Hi 
observations of other galaxies, is a question that should 
be addressed in such a model. 

These new results obtained for M31 could only be ob- 
tained by a combination of deep exposures with high 
spectral resolution measurements. The quest of under- 
standing the formation and evolution of M31 should ben- 
efit from such deep observations. We expect that future 
radio observatories like the Square Kilometre Array (in 
2016) and its australian and south-african precursors (in 
2013) will renew the vision we have of the outer Hi distri- 
bution and dynamics for millions of galaxies in the local 
Universe. For the first time they will also allow to trace 
the evolution of the atomic gas as a function of redshift 
and thus the growth of the baryonic mass of Milky Way- 
and M31-like galaxies. 
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TABLE 4 

Mass model results using the i?— band photometry. 
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TABLE 5 

Hi rotation curve, surface density, disc inclination and position angle from tilted-ring model results. 



Radius 


Radius 


P.A. 


ep A 




i 




7 

''a.doptscl 




^rot 




(') 


Cknc") 


V ) 






V / 


\ } 


\ } 


(km s~l) 


(km s~l) 


{ A^r~\ DC 


1.67 


0.38 






26.6 






31.0 






2.77 


3.33 


0.76 






26.6 






31.0 






3.19 


5.00 


1.14 






26.6 






31.0 


336.2 


171.7 


4.21 


6.67 


1.52 






26.6 






31.0 


324.6 


125.1 


3.94 


8.33 


1.90 


26.6 


2.5 


26.8 


30.5 


4.3 


31.0 


339.0 


52.8 


4.49 


10.00 


2.28 


27.1 


2.2 


27.3 


62.4 


4.3 


51.0 


243.6 


25.8 


2.14 


11.67 


2.66 


27.7 


1.9 


28.6 


61.0 


3.4 


61.3 


235.2 


17.0 


1.52 


13.33 


3.04 


30.7 


1.5 


30.6 


66.3 


2.7 


65.8 


238.9 


5.7 


1.51 


15.00 


3.43 


33.5 


1.3 


32.2 


70.4 


3.0 


68.1 


239.3 


18.3 


1.57 


16.67 


3.81 


32.6 


1.2 


32.7 






67.7 


226.3 


16.1 


2.21 


18.33 


4.19 


32.3 


1.0 


32.9 


69.3 


1.5 


67.3 


202.6 


4.7 


3.00 


20.00 


4.57 






33.3 


62.5 


1.0 


64.4 


207.3 


10.7 


3.54 


21.67 


4.95 






33.6 


60.8 


1.2 


63.7 


202.5 


21.7 


2.75 


23.33 


5.33 


33.4 


0.9 


33.9 


66.8 


1.3 


65.9 


208.9 


15.6 


2.31 


25.00 


5.71 


36.0 


0.5 


35.4 


69.7 


1.4 


68.1 


221.6 


13.4 


2.09 


26.67 


6.09 


37.0 


0.6 


36.4 


66.8 


0.7 


69.7 


232.2 


13.7 


1.86 


28.33 


6.47 


36.3 


0.5 


36.6 


74.3 


0.8 


72.0 


237.6 


8.3 


1.41 


30.00 


6.85 


36.9 


0.4 


36.5 


73.2 


0.4 


73.5 


239.8 


2.2 


1.43 


31.67 


7.23 


35.7 


0.3 


36.4 


74.6 


0.6 


74.3 


235.6 


6.1 


1.60 


33.33 


7.61 


37.0 


0.2 


36.8 


75.0 


0.4 


74.6 


241.7 


3.3 


2.07 


35.00 


7.99 


37.2 


0.2 


37.3 


74.6 


0.4 


74.5 


244.3 


6.4 


2.84 


36.67 


8.37 


37.9 


0.2 


37.7 


73.7 


0.4 


74.3 


248.8 


6.4 


3.85 


38.33 


8.75 


38.1 


0.3 


38.0 


74.8 


0.4 


74.3 


251.8 


5.5 


4.61 


40.00 


9.13 


38.0 


0.3 


38.3 


73.8 


0.5 


74.4 


253.0 


9.2 


4.39 


41.67 


9.51 


38.9 


0.2 


38.7 


75.0 


0.5 


74.8 


258.8 


9.6 


3.78 


43.33 


9.90 


39.1 


0.2 


39.0 


75.2 


0.4 


75.2 


259.0 


9.5 


3.34 


45.00 


10.28 


39.2 


0.2 


39.1 


75.3 


0.4 


75.6 


262.2 


10.8 


3.16 


46.67 


10.66 


39.1 


0.2 


39.0 


76.7 


0.4 


76.1 


266.8 


13.0 


3.37 


48.33 


11.04 


38.7 


0.2 


38.8 


76.1 


0.3 


76.3 


266.8 


11.7 


4.34 


50.00 


11.42 


38.5 


0.2 


38.6 


76.6 


0.3 


76.4 


265.9 


9.9 


5.02 


51.67 


11.80 


38.7 


0.2 


38.3 


76.5 


0.3 


76.3 


264.4 


7.2 


4.95 


53.33 


12.18 


37.8 


0.2 


37.7 


76.1 


0.3 


76.0 


264.7 


5.3 


4.45 


55.00 


12.56 


36.8 


0.2 


37.0 


75.5 


0.3 


75.6 


265.3 


5.4 


4.03 


56.67 


12.94 


36.1 


0.2 


36.4 


75.1 


0.3 


75.1 


265.2 


3.9 


3.98 


58.33 


13.32 


36.3 


0.2 


36.1 


74.7 


0.3 


74.8 


262.0 


2.4 


3.91 


60.00 


13.70 


35.5 


0.2 


35.9 


74.7 


0.3 


74.5 


260.8 


4.5 


3.67 


61.67 


14.08 


36.0 


0.2 


36.1 


73.9 


0.4 


74.3 


259.2 


5.4 


3.39 


63.33 


14.46 


36.6 


0.2 


36.6 


74.9 


0.3 


74.1 


258.1 


6.8 


3.37 


65.00 


14.84 


37.3 


0.2 


37.1 


72.9 


0.3 


73.7 


258.4 


6.1 


3.01 


66.67 


15.23 


37.3 


0.2 


37.3 


73.7 


0.3 


73.7 


259.2 


4.6 


2.61 


68.33 


15.61 


37.4 


0.2 


37.2 


73.9 


0.3 


74.1 


262.7 


4.2 


2.39 


70.00 


15.99 


36.8 


0.2 


37.1 


74.5 


0.3 


74.7 


266.1 


4.0 


2.00 


71.67 


16.37 


37.2 


0.2 


37.3 


75.8 


0.3 


75.4 


270.0 


2.3 


1.92 


73.33 


16.75 


37.9 


0.2 


37.7 


75.7 


0.3 


75.5 


269.8 


1.0 


1.79 


75.00 


17.13 


38.1 


0.2 


37.7 


75.6 


0.3 


75.1 


269.1 


3.9 


2.21 


76.67 


17.51 


37.1 


0.2 


37.5 


73.7 


0.4 


74.3 


268.5 


7.1 


1.91 


78.33 


17.89 


37.5 


0.2 


37.5 


74.2 


0.5 


73.8 


263.0 


15.0 


1.74 


80.00 


18.27 


37.7 


0.3 


37.5 


72.8 


0.5 


73.4 


257.1 


13.5 


1.60 


81.67 


18.65 


37.1 


0.3 


37.4 


73.4 


0.5 


73.4 


254.1 


10.6 


1.44 


83.33 


19.03 






37.5 


73.8 


0.5 


73.6 


251.9 


11.1 


1.35 


85.00 


19.41 


37.1 


0.3 


37.7 


73.7 


0.5 


73.4 


249.5 


8.1 


1.24 


86.67 


19.79 


38.8 


0.3 


38.3 






73.2 


245.7 


7.4 


1.15 


88.33 


20.18 


39.2 


0.3 


38.6 


72.5 


0.6 


73.0 


243.7 


6.6 


1.18 


90.00 


20.56 


37.9 


0.3 


38.3 






73.0 


245.9 


7.5 


1.10 


91.67 


20.94 


37.7 


0.3 


38.2 






73.0 


242.3 


6.1 


1.08 


93.33 


21.32 


38.9 


0.3 


38.4 


73.0 


0.5 


73.0 


239.2 


6.3 


0.83 


95.00 


21.70 


38.3 


0.3 


38.2 


73.0 


0.6 


73.3 


239.5 


4.7 


0.75 


96.67 


22.08 


37.5 


0.3 


37.8 


74.3 


0.6 


73.5 


236.1 


1.8 


0.76 


98.33 


22.46 


37.4 


0.3 


37.5 


72.9 


0.4 


73.4 


233.8 


1.7 


0.88 


100.00 


22.84 


37.5 


0.3 


37.6 


73.3 


0.5 


73.3 


233.1 


3.3 


0.82 


101.67 


23.22 


38.0 


0.3 


37.7 


73.6 


0.5 


73.3 


230.1 


5.5 


0.75 


103.33 


23.60 


37.4 


0.3 


37.7 


72.9 


0.5 


73.2 


232.1 


5.0 


0.67 


105.00 


23.98 


37.6 


0.3 


37.7 


73.4 


0.6 


73.2 


228.7 


1.7 


0.58 


106.67 


24.36 






37.8 


73.0 


0.6 


73.2 


229.1 


1.8 


0.49 


108.33 


24.75 


38.0 


0.3 


38.0 


73.5 


0.6 


73.4 


227.9 


5.0 


0.43 


110.00 


25.13 






38.0 


73.8 


0.6 


73.6 


226.9 


2.0 


0.42 


111.67 


25.51 


38.5 


0.3 


38.1 






73.6 


225.1 


1.6 


0.45 


113.33 


25.89 


37.6 


0.3 


38.1 






73.7 


225.4 


1.8 


0.46 


115.00 


26.27 


38.4 


0.3 


38.3 


73.1 


0.5 


73.7 


230.3 


2.0 


0.42 


116.67 


26.65 






38.4 


74.5 


0.5 


74.0 


229.0 


2.3 


0.39 


118.33 


27.03 


38.9 


0.3 


38.5 


74.2 


0.5 


73.9 


229.9 


4.8 


0.38 


120.00 


27.41 


38.0 


0.3 


38.4 


73.3 


0.7 


73.5 


230.1 


6.6 


0.29 


121.67 


27.79 


38.5 


0.3 


38.2 


72.4 


0.9 


73.6 


229.8 


3.0 


0.23 


123.33 


28.17 


37.8 


0.3 


37.8 


75.6 


0.5 


74.0 


230.4 


5.2 


0.20 


125.00 


28.56 


37.2 


0.3 


37.4 


73.7 


0.9 


73.5 


230.9 


2.9 


0.21 


126.67 


28.94 


37.0 


0.3 


37.0 


71.5 


1.1 


72.5 


229.8 


2.1 


0.17 


128.33 


29.32 


36.6 


0.5 


36.4 


73.3 


1.5 


71.3 


228.8 


1.8 


0.10 


130.00 


29.70 


36.2 


0.4 


34.9 


68.2 


2.0 


69.3 


238.3 


3.3 


0.08 



